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Abstract 

In a recent paper it was shown that, for chemical reaction networks 

. 2 ' possessing a subtle structural property called concordance, dynamical 

,^ ' behavior of a very circumscribed (and largely stable) kind is enforced, 

Q^' so long as the kinetics lies within the very broad and natural weakly 

monotonic class. In particular, multiple equilibria are precluded, as 
ff^ I are degenerate positive equilibria. Moreover, under certain circum- 

^ ' stances, also related to concordance, all real eigenvalues associated 

vQ . with a positive equilibrium are negative. Although concordance of a 

l/^ I reaction network can be decided by readily available computational 

^^ ' means, we show here that, when a nondegenerate network's Species- 

C^ . Reaction Graph satisfies certain mild conditions, concordance and its 

dynamical consequences are ensured. These conditions are weaker 
than earlier ones invoked to establish kinetic system injectivity, which, 
in turn, is just one ramification of network concordance. Because the 
Species- Reaction Graph resembles pathway depictions often drawn by 
biochemists, results here expand the possibility of inferring significant 
C^ ■ dynamical information directly from standard biochemical reaction 

diagrams. 
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1 Introduction 

1.1 Background 

This article is intended as a supplement to another one [18], in which we 
defined the large class of concordant chemical reaction networksi}] and deduced 
many of the striking properties common to all members of that class. We 
argued that, so long as the kinetics is weakly monotonic ( §3.2p . network 
concordance enforces behavior of a very circumscribed kind. 

Among other things, we showed that the class of concordant networks 
coincides precisely with the class of networks which, when taken with any 
weakly monotonic kinetics, invariably give rise to kinetic systems that are 
injective — a quality that precludes, for example, the possibility of switch- 
like transitions between distinct positive stoichiometrically compatible steady 
states. Moreover, we showed that certain properties of concordant networks 
taken with weakly monotonic kinetics extend to still broader categories of 
kinetics — including kinetics that involve product inhibition — provided that 
the networks considered are not only concordant but also strongly concordant 



Although reaction network concordance is a subtle structural property, 
determination of whether or not a network is concordant (or strongly concor- 
dant) is readily accomplished with the help of easy-to-use and freely available 
software [15], at least if the network is of moderate sizej3 In this way, one can 
easily determine whether the several dynamical consequences of concordance 
or strong concordance accrue to a particular reaction network of interest. 

1.2 The role of the Species-Reaction Graph 

The Species-Reaction Graph (SR Graph), defined in Section [2|, is a graphical 
depiction of a chemical reaction network resembling pathway diagrams often 
drawn by biochemists. For so-called "fully open" systems ( §1.2.ip . earlier 
work [Il|2l[71l9] indicated that, when the SR Graph satisfies certain structural 
conditions and when the kinetics is within a specified class, the governing 
differential equations can only admit behavior of a restricted kind. A survey 
of some of that earlier work, beginning with the Ph.D. research of Paul 



^Thc formal definitions of concordance and strong concordance appear here in Sections 
Hand El 

^Some of the underlying algorithms are given in |14j . 



Schlosser [inillZ] is provided in [IB] . Although most of the initial SR Graph 
results for fully open networks were focused on mass action kinetics, that 
changed with the surprising work of Banaji and Craciun [Ill2]- (In contrast 
to the SR Graph results, however, the network analysis tools in [15] are 
indifferent to whether the network is fully open.) 

With this as background, in pES] we asserted without proof that attributes 
of the SR Graph shown earlier to be sufficient for other network properties 
[Ill2l[71[9] are, in fact, sufficient to ensure not only concordance but also 
strong concordance, at least in the "fully open" setting. In turn, those other 
properties (e.g., the absence of multiple steady states when the kinetics is 
weakly monotonic) largely derive from concordance. That certain SR Graph 
attributes might ensure concordance seems, then, to be the fundamental idea, 
with attendant dynamical consequences of those same SR Graph attributes 
ultimately descending from concordance. 

Although handy computational tools in [15] will generally be more inci- 
sive than the SR Graph in determining a network's concordance properties 
and although those computational tools are indifferent to whether the net- 
work is fully open, the SR Graph nevertheless has its strong attractions. 
Not least of these is the close relationship that the SR Graph bears to re- 
action network diagrams often drawn by biochemists. For this reason alone, 
far-from-obvious dynamical consequences that might be inferred from these 
ubiquitous diagrams take on considerable interest. Moreover, theorems that 
tie network concordance to properties of the SR Graph point to consequences 
of reaction network structure that are not easily gleaned from implementation 
of computational tests on a case-by-case basis. 

Thus, our primary purpose here is to prove the assertions made in [18] 
about connections between concordance of a network and the nature of its 
SR Graph. In fact, we go further in two directions, directions that would 
have been somewhat askew to the main thrusts of 1181: 



1.2.1 The Species-Reaction Graph and the concordance of "fully 
open" networks 

First, we show that one can infer concordance of a so-called "fully open" 
network when its SR graph satisfies conditions, stated in Theorem 12. 1^ that 
are substantially weaker than those stated in [TSj (or in [Il[21[71in] ) ; Corollary 
12.41 invokes the SR Graph conditions stated in [18]. A fully open network is 
one containing a "degradation reaction" s -> for each species s in the net- 



workj^ Proof of concordance when the aforementioned SR Graph conditions 
are satisfied is by far the main undertaking of this article. 

1.2.2 The Species-Reaction Graph and the concordance of net- 
works that are not "fully open" 

Second, we examine more completely concordance information that the SR 
Graph gives for reaction networks that are not fully open. 

For a network that is not fully open, its fully open extension is the net- 
work obtained from the original one by adding degradation reactions for all 
species for which such degradation reactions are not already present. In [18] 
we proved that a normal network is concordant (respectively, strongly con- 
cordant) if its fully open extension is concordant (strongly concordant). A 
normal network is one satisfying a very mild condition described in [8], [T8] . 
and, here again, in §7.11 The large class of normal networks includes all 
weakly reversible networks and, in particular, all reversible networks /^. 

Thus, if inspection of the SR Graph can serve to establish, by means of 
theorems crafted for fully open networks, that a normal network's fully open 
extension is concordant or strongly concordant, then the original network 
itself has that same property. In this way, the seemingly restricted power of 
"fully open" SR graph theorems extends far beyond the fully open setting. 

This has some importance, especially in relation to reversible networks. 
Consider a reaction network whose SR Graph ensures, by virtue of Corollary 
12. 4[ concordance of the network's fully open extension. Then the original 
network can fail to be concordant only if it is not reversible. But a kinetic 
system based on a network that is not reversible is usually deemed to be an 
approximation of a more exact "nearby" kinetic system in which all reactions 
are reversible, perhaps with some reverse reactions having extremely small 
rates. 

To the extent that this is the case, lack of concordance in the original 
network is an artifact of the approximation: The reversible network under- 
lying the more exact kinetic system is concordant, so that system inherits 
all of the dynamical properties in [T^ that are consequences of network con- 
cordance. (In particular, it inherits those properties listed as items (i) - (iii) 
in Theorem 12.81 ) Thus, any absence of these attributes in a kinetic system 
based on the original network is, again, an artifact of the approximation. 



■^Nevertheless, the SR graph is drawn for the so-called "true" chemistry, devoid of 
degradation-synthesis reactions such as s ^ or ^ s. 
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Even so, kinetic system models based on networks that are not reversible 
(or, more generally, weakly reversible) have an intrinsic interest. They are 
ubiquitous, and one would like to understand their inherent properties. Re- 
sults in y^ about network normality and concordance already go a long way 
in this direction, for the class of normal networks extends well beyond the 
weakly reversible class. 

Here we argue that these same results extend to the still larger class of 
nondegenerate reaction networks, a class that includes all normal networks 
and, in particular, all weakly reversible networks. A network is nondegener- 
ate if, taken with some differentiably monotonic kinetics [18], the resulting 
kinetic system admits even one positive composition at which the derivative 
of the species-formation-function is nonsingulaio; otherwise we say that the 
network is degenerate (Definition 17. 6p . 

In fact, degenerate reaction networks are never concordant ( §7.ip . so for 
them questions about the possibility of concordance are moot. On the other 
hand, a nondegenerate network is concordant (strongly concordant) if the 
network's fully open extension is concordant (strongly concordant) ( §7.2p . 
Thus, SR Graph theorems that give information about the concordance or 
strong concordance of a nondegenerate network's fully open extension give 
that same information about the network itself. 

Implications of network normality were considered in some depth in [18] . 
Because corresponding considerations of the broader notion of nondegeneracy 
are similar and because those considerations are rather different in spirit 
from the largely graph-theoretical aspects of most of this article, we chose 
to defer the entire discussion of nondegeneracy to Section [71 which includes 
computational tests whereby normality and nondegeneracy of a network can 
be affirmed. 

In any case, it should be kept in mind that the computational tools pro- 
vided in [15] are indifferent to whether the network under study is fully open 
(or nondegenerate). 



^A positive composition is one at which all species concentrations are (strictly) positive. 
When we say that the derivative ( Jacobian matrix) of the species-formation rate function 
is nonsingular we mean that its null space contains no nonzero vector of the stoichiometric 
subspace. 



1.3 The Species- Reaction Graph and consequences of 
concordance 

For the most part the theorems in this article have as their objective the 
drawing of connections between the concordance of a reaction network and 
properties of the network's Species- Reaction Graph. However, it is important 
to keep in mind the dynamical consequences of these theorems. When an 
SR Graph theorem asserts concordance of a particular network, the same 
theorem is also asserting that the network inherits all of the properties shown 
in fT8^ to accrue to all concordant networks. 

Thus, for example, a theorem that asserts concordance for any reaction 
network whose SR Graph has Property X is also asserting that, for any 
reaction network whose SR Graph has Property X, there is no possibility 
of two distinct stoichiometrically compatible equilibria, at least one of which 
is positive, no matter what the kinetics might be, so long as it is weakly 
monotonic. 

A theorem that does make explicit connections between the SR Graph 
and dynamical consequences is offered at the close of the next section. 

2 Main theorems 

The Species- Reaction Graph (SR Graph) of a chemical reaction network is 
a bipartite graph constructed in the following way: The vertices are of two 
kinds — species vertices and reaction vertices. The species vertices are simply 
the species of the network. The reaction vertices are the reactions of the 
network but with the understanding that a reversible reaction pair such as 
A + B :^ P is identified with a single vertex. If a species appears in a 
particular reaction, then an edge is drawn that connects the species with 
that reaction, and the edge is labeled with the name of the complex in which 
the species appears. (The complexes of a reaction network are the objects 
that appear before and after the reaction arrows. Thus, the complexes of 
reaction A + B ^ P are A + B and P.) The stoichiometric coefficient of an 
edge is the stoichiometric coefficient of the adjacent species in the labeling 
complex. We show the Species- Reaction Graph for network ([T]) in Figure [TJ 
The arrows on some of the edges will be explained shortly. 
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Figure 1: An example of a Species- Reaction Graph 

It will be understood that our focus is on concordance of a "fully open" 
reaction network {=5^,"^,^}. To say that the network is fully open is to 
say that, for each species s in the network there is a "degradation reaction" 
of the form s -> 0. (There might also be "synthesis reactions" of the form 
^ s.) Moreover, we shall also suppose hereafter that a species can appear 
on only one side of a reaction. The first of these requirements was invoked 
in [21E1E], while the second was invoked only in [2]. 



Even in consideration of fully open reaction networks we will always re- 
strict our attention to the Species-Reaction Graph for the "true chemistry" 
— that is, the Species- Reaction Graph for the original fully open network but 
with degradation-synthesis "reactions" such as A ^ or ^ v4 omitted. As 
a reminder we will sometimes refer to the "true-chemistry Species-Reaction 
Graph," but this understanding will always be implicit even when it is not 
made explicit. 

By the intersection of two cycles in the Species- Reaction Graph, we mean 
the subgraph consisting of all vertices and edges common to the two cycles. 
We say that two cycles have a species-to-reaction intersection (S-to-R inter- 
section) if their intersection is not empty and each of the connected compo- 
nents of the intersection is a path having a species at one end and a reaction 
at the other. In the phrase "S-to-R intersection" no directionality is implied. 
In Figure [H Cycle I and Cycle II have an S-to-R intersection consisting of the 
edge connecting species C to reactions B + C <^ Q. There is a third unlabeled 
"outer" cycle, hereafter called Cycle III, that traverses species A,C,Q,B and 
returns to A. It has an S-to-R intersection with Cycle I, in particular, the 
long outer path connecting reactions i? + C «± Q to species C via species A. 
Cycle III also has an S-to-R intersection with Cycle II. 

An edge-pair in the SR Graph is a pair of edges adjacent to a common 
reaction vertex. A c-pair (abbreviation for complex-pair) is an edge-pair 
whose edges carry the same complex label. For readers with access to color, 
the c-pairs in Figure [T] are given distinct colors. 

An even cycle in the SR Graph is a cycle whose edges contain an even 
number of c-pairsu Cycle I contains two c-pairs and Cycle II contains none, 
so both are even. Cycle III, the large outer cycle, contains one c-pair, so it 
is not even. 

A fixed- direction edge-pair is an edge-pair that is not a c-pair and for 
which the reaction of the edge-pair is irreversible. In this case, we assign 
a fixed direction to each of the two edges in the intuitively obvious way: 
The edge adjacent to the reactant species (i.e., the species appearing in the 
reactant complex) is directed away from that species, and the edge adjacent 
to the product species (i.e., the species appearing in the product complex) 
points toward that species. Thus, for example, a fixed-direction edge-pair 



^In the graph theory hterature, the term "even cycle" sometimes refers to a cycle 
containing an even number of edges. That is never the usage here. In a bipartite graph 
— in particular in an SR graph — cycles always have an even number of edges. 
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In Figure [T] there are fixed-direction edge-pairs centered at the reactions 
2A -^ C and C + D -^ Q + E; the corresponding fixed-directions for the 
adjacent edges are shown in the figure. 

An orientation for a simple cycle siRiS2R2---SnRnSi in the SR Graph is 
an assignment of one of two directions to the edges, either 

Si -> i?l ^ S2 -^ R2--- ^ Sn-^ Rn ^ «! 

or 

Si -^ Rl ^ S2 ^ R2--- ^ Sn^ Rn ^ Si, 

that is consistent with any fixed directions the cycle might contain. A cycle 
is orientable if it admits at least one orientation. By an oriented cycle in the 
SR Graph we mean an orientable cycle taken with a choice of orientation. 

Note that a cycle that has no fixed-direction edge-pair will be orientable 
and have two orientations. (This happens when every reaction in the cycle 
is either reversible or else is contained in a c-pair within the cycle.) A cycle 
that has just one fixed-direction edge-pair is orientable and has a unique 
orientation. A cycle that has more than one fixed-direction edge-pair might 
or might not be orientable, depending on whether the fixed-direction edge- 
pairs all point "clockwise" or "counterclockwise," but if there is an orientation 
it will be unique. All cycles in Figured] are orientable, with each cycle having 
a unique orientation. 

A set of cycles in an SR graph admits a consistent orientation if the cycles 
can all be assigned orientations such that each edge contained in one of the 
cycles has the same orientation-direction in every cycle of the set in which 
that edge appears. To see that a set of cycles might not admit a consistent 
orientation, consider a pair of cycles sharing a common edge, each cycle 
having the same compulsory orientation (e.g., "clockwise"). In fact, this is 
precisely the situation for Cycles I and II, so the set consisting of those two 
cycles does not admit a consistent orientation. On the other hand the set 
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consisting of Cycles I and III (the large outer cycle) does admit a consistent 
orientation, as does the set consisting of Cycles II and III. 

A critical subgraph of an SR Graph is a union of a set of even cycles taken 
from the SR Graph that admits a consistent orientation. Because Cycle III 
is not even, it is easy to see that there are only two critical subgraphs in 
Figure [TJ One consists only of Cycle I and the other consists only of Cycle 
II. 

Consider an oriented cycle in the SR Graph, si -^ i?i -> S2 -^ -R2 •••-*■ s^ -> 
Rn -^ si. For each directed reaction-to-species edge R -^ s we denote by Jr^s 
the stoichiometric coefficient associated with that edge. For each directed 
species-to-reaction edge, we denote by e^^/j the stoichiometric coefficient 
associated with it. The cycle is stoichiometrically expansive relative to the 
given orientation if 

>1. (2) 



JRl->-S2JR2^S3 ■ ■ ■ JRn^Si 



^Si->-Ri(^S2->-R2 ■ ■ ■ ^S„->-Rn 

One of the major goals of this article is proof of the following theorem: 

Theorem 2.1. A fully open reaction network is concordant if its true- chemistry 
Species-Reaction Graph has the following properties: 

(i) No even cycle admits a stoichiometrically expansive orientation. 

(a) In no critical subgraph do two even cycles have a species-to-reaction 
intersection. 

Remark 2.2. If there are no critical subgraphs — in particular if there are 
no orientable even cycles — then the conditions of the theorem are satisfied 
trivially, and the network is concordant. 

Remark 2.3. Determination of which subgraphs of the Species- Reaction Graph 
are critical requires consideration of orientability. However, it should be un- 
derstood that condition (m) of Theorem 12.11 is imposed on every undirected 
critical subgraph of the undirected Species-Reaction Graph. In particular, 
account should be taken of all even cycles in the subgraph, whether or not 
they be orientable. (Although, a critical subgraph is the union of a set of 
even cycles that admit a consistent orientation, the resulting subgraph might 
also contain cycles that are not orientable.) Again, no directionality should 
be associated with the phrase "species-to-reaction" intersection. 

With regard to condition (i), note that a cycle that admits two orienta- 
tions can be non-expansive relative to both orientations only if, relative to 
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one of the orientations, the number on the left side of (|2]) is equal to one 
(in which case it will also have that same value with respect to the opposite 
orientation). 

In the language of [2l[6l[9] , we say that a (not necessarily orientable) cycle 
in the SR Graph is an s-cycle if relative to an arbitrarily imposed "clockwise" 
direction, the calculation on the left side of ([2]) yields the value one. With 
this in mind we can state an "orientation-free" corollary of Theorem 12.11 

Corollary 2.4. A fully open reaction network is concordant if its true- 
chemistry Species- Reaction Graph has the following properties: 

(i) Every even cycle is an s-cycle. 

(a) No two even cycles have a species-to-reaction intersection. 

Especially for networks in which there are several irreversible reactions. 
Theorem 12.11 is likely to be more incisive than its corollary. Consider, for 
example, the Species- Reaction Graph for network ([1]). Both conditions of 
Corollary 12.41 fail. Nevertheless, both oi the (weaker) conditions of Theorem 
12.11 are satisfied, so a fully open network having ([1]) as its true chemistry is 
concordant. 

Remark 2.5. We shall see very easily in Section [7] that network ([T]) is normal. 
Because every normal network with a concordant fully open extension is itself 
concordant, we can conclude that network ([1]), like its fully open extension, 
is concordant. 

When a fully open network does satisfy the stronger conditions of Corol- 
lary |2]1] we can say more, this time about strong concordance. 

Theorem 2.6. A fully open reaction network is strongly concordant if its 
true- chemistry Species- Reaction Graph has the following properties: 

(i) Every even cycle is an s-cycle. 

(a) No two even cycles have a species-to-reaction intersection. 

Although the fully open extension of network ([1]) is concordant, it is not 
strongly concordant. (This can be quickly determined by the freely available 
Ghemical Reaction Network Toolbox [15j.) The Species- Reaction Graph for 
network ([T]) does not satisfy the conditions of Theorem | 
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Remark 2.7. Although these theorems nominally speak in terms of concor- 
dance or strong concordance of a fully open network, they are even more 
generous than they seem: Recall that in Section [7] we indicate why the con- 
cordance properties ensured by the theorems actually extend beyond the fully 
open setting to a far wider class of reaction networks. In particular, they 
extend to the class of nondegenerate networks described earlier in §1.2.21 (in- 
cluding all weakly reversible networks). For nondegenerate networks, then, 
the fully open requirement of Theorems 12.11 and 12.61 becomes moot. 

In the spirit of §1.31 we close this section with the statement of a theo- 
rem that ties the hypothesis of Theorem 12.11 not to concordance itself but, 
instead, to dynamical consequences of concordance proved in |18]. The fol- 



lowing theorem is essentially a corollary of Theorem 12. 6[ theorems in [18] . 
and material to be discussed in Section [71 Language in the theorem is that 
used in ^18j. Much of it is reviewed in Section E] of this article. 

Theorem 2.8. For any nondegenerate reaction network whose true chem- 
istry Species-Reaction Graph satisfies conditions (i) and (ii) of Theorem \2.1\ 
the following statements hold true: 

(i) For each choice of weakly monotonic kinetics, every positive equilibrium is 
unique within its stoichiometric compatibility class. That is, no positive equi- 
librium is stoichiometrically compatible with a different equilibrium, positive 
or otherwise. 

(ii) If the network is weakly reversible then, for each choice of kinetics (not 
necessarily weakly monotonic), no nontrivial stoichiometric compatibility class 
has an equilibrium on its boundary. In fact, at each boundary composition in 
any non-trivial stoichiometric compatibility class the species-formation-rate 
vector points into the stoichiometric compatibility class in the sense that there 
is an absent species produced at strictly positive rate. If, in addition, the net- 
work is conservative then, for any choice of a continuous weakly monotonic 
kinetics, there is precisely one equilibrium in each nontrivial stoichiometric 
compatibility class, and it is positive. 

(Hi) If the kinetics is differentiably monotonic then every positive equilibrium 
is nondegenerate. Moreover, every real eigenvalue associated with a positive 
equilibrium is negative. 
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As in [TS], it is understood that eigenvalues in the theorem statement are 
those associated with eigenvectors in the stoichiometric subspace. Similarly, 
when we say that a positive equilibrium in nondegenerate, we mean that, for 
the equilibrium, zero is not an eigenvalue corresponding to an eigenvector in 
the stoichiometric subspace. A stoichiometric compatibility class is nontrivial 
if it contains at least one positive composition. 

Remark 2.9. Consider a (not necessarily nondegenerate) reaction network 
whose true chemistry SR Graph satisfies conditions (i) and (ii) of Theorem 
12. 1[ With respect to the possibility of nondegenerate positive equilibria. The- 
orem 12.81 describes something of an all or nothing situation: If there is some 
different iably monotonic kinetics that gives rise to even one nondegenerate 
positive equilibrium, then the network itself is nondegenerate, in which case 
every positive equilibrium arising from any differentiably monotonic kinetics 
is nondegenerate (and unique within its stoichiometric compatibility class). 

Remark 2.10. Consider a nondegenerate reaction network that is not neces- 
sarily fully open. If its true-chemistry SR Graph satisfies the conditions of 
Theorem 12. 6[ then considerations in Section [7] will indicate that the network 
is strongly concordant. In this case, the dynamical properties of strongly con- 
cordant networks given in [18] accrue to the network at hand, and one can 
again deduce a theorem which, in the spirit of Theorem 12. 8 [ makes statements 
about general properties of kinetic systems the network engenders, this time 
including those that derive from so-called "two-way monotonic kinetics." 

3 Reaction network theory preliminaries 

This section is, for the most part, is a compendium of the definitions and 
infrastructure used in [18], repeated here for the reader's convenience. Ref- 
erence [18j has more in the way of discussion, and [10] has still more in the 
way of motivation. We begin with notation. 

3.1 Notation 

When J is a finite set (for example, a set of species or a set of reactions), we 
denote the vector space of real- valued functions with domain J by M-^. If A^ 
is the number of elements in the set /, then W is, in effect, a copy of the 
standard vector space M^, with the components of a vector x e M^ indexed by 
the names of the members of / instead of the integers 1, 2, . . . , A. For x € M-^ 

15 



and i e I, the symbol Xi denotes the value (component) of x corresponding 
to the element i e I. 

For example, if ^ = {NO,02,N02} is the set of species in a chemical 
system and if cno,co2, and C7VO2 are the molar concentrations of the three 
species in a particular mixture state, then that state can be represented by 
a "composition vector" c in the vector space M.-^ . That is, c represents an 
assignment to each species of a number, the corresponding molar concentra- 
tion. Readers who wish to do so can, in this case, simply regard c to be the 
3- vector [cNo,co2,CN02]y with the understanding that the arrangement of the 
three numbers in such an ordered array is superfluous to the mathematics at 
hand. 

Vector representations in spaces such as W rather than M-^ have advan- 
tages in consideration of graphs and networks, where one wants to avoid 
nomenclature that imparts an artificial numerical order to vertices or edges. 
See, for example, |3]. 

The subset of W consisting of vectors having only positive (nonnegative) 

components is denoted Rl (1R+)- For each x € M-^ and for each 2; e M^, the 
symbol x^ denotes the real number defined by 



x':-- 






where it is understood that 0" = 1. For each x,x' e W, the symbol x o x' 
denotes the vector of MJ defined by 

(x o x')i ■■= XiX'i, Vi € I. 

For each i € /, we denote by Ui the vector of W such that {u!i)j = 1 whenever 
j = i and (u}i)j = whenever j i^ i. 

The standard basis for W is the set {ui € W : i e 1} . Thus, for each x € W , 
we can write x = Y,iei^i^i- The standard scalar product in W is defined as 
follows: If X and x' are elements of M^, then 

tjO Ju — / ih '^Ju A t 

The standard basis of M^ is orthonormal with respect to the standard scalar 
product. It will be understood that MJ carries the standard scalar product 
and the norm derived from the standard scalar product. It will also be 
understood that M^ carries the corresponding norm topology. 
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If f/ is a linear subspace of R^, we denote by [/-"- the orthogonal comple- 
ment of U in W with respect to the standard scalar product. 

By the support of x € M.^ , denoted suppa;, we mean the set of indices z € / 
for which Xi is different from zero. When ^ is a real number, the symbol 
sgn (^) denotes the sign of ^. When x is a vector of W , sgn (x) denotes the 
function with domain / defined by 

(sgn(x))i :=sgna;i, Vie/. 

3.2 Some definitions 

As indicated in Section [2], the objects in a reaction network that appear at 
the heads and tails of the reaction arrows are the complexes of the network. 
Thus, in network ([2]), the complexes are 2A,B,C,C + D and E. A reaction 
network can then be viewed as a directed graph, with complexes playing the 
role of the vertices and reaction arrows playing the role of the edges. 

2A ^ B 

C (3) 

C + D^E 

Remark 3.1. Let ^ be the set of species in a network. In chemical reaction 
network theory, it is sometimes the custom to replace symbols for the stan- 
dard basis of R'^ with the names of the species themselves. For example, in 
network ([3]), with y = {A,B,C,D,E}, a vector such as uc + ojd e R'^ can 
instead be written as C + D, and 2ua can be written as 2A. In this way, R-^ 
can be identified with the vector space of formal linear combinations of the 
species. As a result, the complexes of a reaction network with species set 
y can be identified with vectors in R'^. 

Definition 3.2. A chemical reaction network consists of three finite sets: 

1. a set =5^ of distinct species of the network; 

2. a set ^ c R:^ of distinct complexes of the network; 

3. a set .^ c "^ X "^ of distinct reactions, with the following properties: 
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(a) {y, y)iM for any y e ^; 

(b) for each y ^'^ there exists y' e ^ such that {y,y') e ^ or such 
that (y',y) e ^. 

If {y,y') is a member of the reaction set ^, we say that y reacts to y', 
and we write y -^ y' to indicate the reaction whereby complex y reacts to 
complex y'. The complex situated at the tail of a reaction arrow is the 
reactant complex of the corresponding reaction, and the complex situated at 
the head is the reaction's product complex. 

For network ©, ^ = {A,B,C,D,E}, "€ = {2A,B,C,C + D,E}, and 
^ = {2A ^ B,B ^ 2A,B ^ C,C ^ 2A,C + D ^ E,E ^ C + D}. The 
diagram in (^ is an example of a standard reaction diagram: each complex 
in the network is displayed precisely once, and each reaction in the network 
is indicated by an arrow in the obvious way. 

In the context of the present paper and its predecessor ^Bj the idea of 
weak reversibility plays an important role. The following definition provides 
some preparation. 

Definition 3.3. A complex y ^ ^ ultimately reacts to a complex y' e ^ if 
any of the following conditions is satisfied: 

1. y^y'e^g; 

2. There is a sequence of complexes y{l),y{2), . . . ,y(k) such that 

y^y{l)^y{2)^...^y{k)^y'. 

In network ([3]) the complex 2A ultimately reacts to the complex C, but 
the complex C does not ultimately react to the complex C + D. 

Definition 3.4. A reaction network {=5^,^,^} is reversible ii y' -^ y e ^ 
whenever y ^ y' ^ !!%. The network is weakly reversible if for each y,y' € ^, 
y' ultimately reacts to y whenever y ultimately reacts to y' . 

Network ([3]) is weakly reversible but not reversible. On the other hand, 
every reversible network is also weakly reversible. A reaction network is 
weakly reversible if and only if in its standard reaction diagram every arrow 
resides in a directed arrow-cycle. 



Definition 3.5. The reaction vectors for a reaction network {S^,'io,^} are 
the members of the set 

{y'-yem'^':y^y'e3?}. 

The rank of a reaction network is the rank of its set of reaction vectors. 

For network ([3]) the reaction vector corresponding to the reaction 2A -> B 
is B - 2A. The reaction vector corresponding to the reaction C + D ^ E is 
E-C-D. 

Definition 3.6. The stoichiometric suhspace S oia. reaction network {^, ^, ^} 
is the linear subspace of M'^ defined by 

S := span [y' -y^W-^ -.y^y' ^M) . (4) 

The dimension of the stoichiometric subspace is identical to the network's 
rank. The stoichiometric subspace is a proper linear subspace of M:^ when- 
ever the network is conservative: 

Definition 3.7. A reaction network {^,^,^} is conservative if the or- 
thogonal complement S^ of the stoichiometric subspace S contains a strictly 
positive member of M-^': 

S^ n M:^' ^ 0. 

Network ([3]) is conservative: The strictly positive vector [A + 2B + 2C + 
D + 3E) e Wf is orthogonal to each of the reaction vectors of ([3]). 

For a reaction network {^, ^, ^} a mixture state is generally represented 

— .y 
by a composition vector c e M^ , where, for each s € S^, we understand c^ 

to be the molar concentration of s. By a positive composition we mean a 

strictly positive composition — that is, a composition in Mf^. 

Definition 3.8. A kinetics Jt^ for a reaction network {^,^,M} is an as- 
signment to each reaction y ->■ y' e ^ oi a rate function ^y^y' '■ '^+ ^ M+ 
such that 

J(fy^y'{c) > if and only if supp?y c suppc. 

Definition 3.9. A kinetic system {S^, ^, M, J(f} is a reaction network {S^ ., ^, M} 
taken with a kinetics J^ for the network. 
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Many of the dynamical consequences of network concordance require that 
the kinetics be weakly monotonic or diff'erentiably monotonic. Both are for- 
mally defined below: 

Definition 3.10. A kinetics J(f for reaction network {^5^,^,^^} is weakly 
monotonic if, for each pair of compositions c* and c**, the following impli- 
cations hold for each reaction y ^ y' e M such that suppy c suppc* and 
suppy c suppc**: 

(i) ,^^y'(c**) > J^y^yi{c*) => fheiQ Is a species s e supply with c** > c*. 

(ii) ^^y'{c**) = J^y^y'(c*) =^ c* * = c* for all s e suppy or else there 
are species s, s' e supp y with c** > c* and c*,* < c*,. 

We say that the kinetic system {^,^,^, J^} is weakly monotonic when its 
kinetics Jt^ is weakly monotonic. 

Definition 3.11. A kinetics J^ for a reaction network {S^,'^,M} is diff'er- 
entiably monotonic at c* e Mf^ if, for every reaction y ^ y' e M, ^y^y'{-) is 
differentiable at c* and, moreover, for each species s e ^, 

^#^(c^)>0, (5) 

dCs 

with inequality holding if and only if s e supp y. A diff'erentiably monotonic 
kinetics is one that is differentiably monotonic at every positive composition. 

When a kinetics J^f for a reaction network {^,'^,^}is differentiably 
monotonic at c* e Mf^, we denote by the symbol V^^j;'(c*) the member of 
M'^ defined by 

V^,.,,(c*) := ^ ^§^(c*).. 

Note that every mass action kinetics is both weakly monotonic and dif- 
ferentiably monotonic. Recall that a mass action kinetics is a kinetics in 
which the rate function corresponding to each reaction y -^ y' takes the form 

J^y^y'{c) = ky^y'cy,yc € M^ , where ky^y' is a positive rate constant for the 
reaction y -^ y'. 
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Definition 3.12. The species formation rate function for a kinetic system 
{y, ^,M, J^} witli stoichiometric subspace 5" is the map / : M+ -^ S defined 

by 

/(C) = Y. ^y^y'{c)W-y). (6) 



Definition 3.13. The differential equation for a kinetic system with species 
formation rate function /(•) is given by 

c = /(c). (7) 

From equations (|1]), (|6]), and ([Tj), it follows that, for a kinetic system 
l^,"^,,^,^}, the vector c will invariably lie in the stoichiometric subspace 

y y 

S for the network {^,^,,^} . Thus, two compositions c € M_^ and d e R_^ 
can lie along the same solution of ([7]) only if their difference c' - c lies in S. 
This motivates the following definition: 

Definition 3.14. Let {^,^,^} be a reaction network with stoichiometric 

y 

subspace S. Two compositions c and c' in M^ are stoichiometrically compat- 
ible if c' - c lies in S. 

For a network {=5^, ^, ^} , the stoichiometric compatibility relation serves 
— y 
to partition R_^ into equivalence classes called the stoichiometric compati- 
bility classes for the network. Thus, the stoichiometric compatibility class 

y 

containing an arbitrary composition c, denoted (c + S") n ]R_|_ , is given by 

(c + >S)nlf = |c'€lf :c'-c€5|. (8) 



y 

The notation is intended to suggest that {c + S) n M^ is the intersection of 

y 

]R_i_ with the parallel of S containing c. A stoichiometric compatibility class 
is nontrivial if it contains at least one (strictly) positive composition. 

Definition 3.15. An equilibrium of a kinetic system {^, ^, M, J^} is a com- 

y 

position c e M+ for which /(c) = 0. A positive equilibrium is an equilibrium 
that lies in Rf^. 



21 



Because compositions along solutions of (j7]) are stoichiometrically com- 
patible, one is typically interested in changes in values of the species for- 
mation rate function that result from changes in composition that are stoi- 
chiometrically compatible. In particular, if /(c*) is the value of the species 
formation rate function at composition c*, then one might be interested in 
the value of /(c) for a composition c very close to c* and stoichiometrically 
compatible with it. Thus, for a kinetic system {^,^,^, ^} with stoichio- 
metric subspace S", with smooth reaction rate functions, and with species 
formation rate function / : M^ ^ S", we will want to work with the derivative 
df{c*) ■■ S ->- S, defined by 

df(c* + 9(j) 



df{c*)a 



dO 



, ^aeS. (9) 

=0 



We say that c* £ Mf^ is a nondegenerate equilibrium if c* is an equilibrium 
and if, moreover, df{c*) is nonsingular. An eigenvalue associated with a 
positive equilibrium c* e Mf^ is an eigenvalue of the derivative df{c*). 

4 Concordant reaction networks 

Here we recall the definition of reaction network concordance [18] . We con- 
sider a reaction network {o$^,^,^} with stoichiometric subspace S c M-^, 
and we let L : M* -^ S" be the linear map defined by 

La= Y, ay^y'{y'-y). (10) 

Definition 4.1. The reaction network {^,^,^}is concordant if there do 
not exist an a € ker L and a nonzero a ^ S having the following properties: 

(i) For each y ^ y' such that cty^y' ^ 0, supp y contains a species s for which 
sgn cr^ ^sgn ay^y'. 

(ii) For each y -^ y' such that cty^y' = 0, ag = for all s e supp?/ or else 
suppy contains species sand s' for which sgn cr^ = -sgn as', both not zero. 

A network that is not concordant is discordant. 

Note that for a fully open network with species set S^ the stoichiometric 
subspace coincides with R-^. The following lemma will prove useful later on: 
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Lemma 4.2. // a fully open network is discordant, it is always possible to 
choose from each reversible pair of non- degradation reactions at least one 
(and sometimes both) of the reactions for removal such that the resulting 
fully open subnetwork is again discordant. 

Proof. Suppose that a fully open network {e5^,^,^}is discordant. Then 
there are a € kerL and nonzero a that together satisfy conditions (i) and 
(n) in Definition 14.11 In particular, we have 

^ ay^yr{y'-y) = 0. (11) 

Let ?/ 4^ y' be a pair of reversible non-degradation reactions in ^, and 
suppose that ay^yf 4^ Ciy'^y with the complexes labeled such that |ag^y'| > 
la^'^yl. In this case, let ^'^ ■■= M \ {y' -^ y}, let a'--, ■■= Oiy^yi - a-g'^y, and, 
for all other y ^ y' in ^\ let cCy^yi '■= oiy^yi. From (TTT]) it follows easily that 

E ^U'W-y) = ^- (12) 

If, on the other hand, ay^yi = cxyi^y, then we can choose ^t := ^\{y' ^^ y}, 
and, for all y -^ y' in ^\ we can let cty , ■■= oty^y'- In this case, f lT2|) will 
again obtain. 

In either case, it is easy to see that a^, taken with the original a, suffices 
to establish the discordance of the subnetwork associated with the reaction 

set ^t. n 

Remark 4.3. In effect. Lemma [4.21 tells us that every fully open discordant 
network with reversible non- degradation reactions has a discordant fully open 
subnetwork in which no non-degradation reaction is reversible. Note that, 
apart from minor changes in labels within the reaction nodes (i.e., replace- 
ment of ?/ ;± I/' by y -^ y') the Species-Reaction Graph derived from the in- 
dicated discordant subnetwork is a subgraph of the Species-Reaction Graph 
drawn for the original network. As we explain at the beginning of Section 
ISl that subgraph will satisfy the conditions of Theorem 12.11 and its corollary 
if the parent Species- Reaction Graph does. These observations will help us 
simplify certain arguments that are otherwise complicated by the presence 
of reversible reactions. 
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5 Proof of Theorem 12.1 



The proof will be by contradiction. With this in mind, we assume hereafter 
the true-chemistry Species-Reaction Graph (SR Graph) for the fully open 
network under consideration has both attributes of the theorem statement 
and that, contrary to the assertion of the theorem, the fully open network is 
discordant. 

In this case. Lemma 14.21 tells us that, when the true chemistry has re- 
versible reaction pairs, then each such pair can be replaced by a certain 
irreversible reaction of the pair (or else removed completely) such that the 
resulting fully open network is again discordant. The SR Graph for the al- 
tered (totally irreversible) true chemistry might have fewer cycles than in the 
original SR Graph (but never more) because of removal of reversible reac- 
tion pairs. Similarly, there might be fewer orientable cycles (but never more) 
than in the original SR Graph as a result of replacement of reversible reaction 
pairs by single irreversible reactions. Moreover, there might be fewer critical 
subgraphs (but never more) than in the original SR Graph because of a loss 
of orientable cycles or because more cycles have acquired compulsory orien- 
tations. (Cycles that were even in the original SR Graph and persist in the 
new SR Graph remain even.) For all of these reasons, the SR Graph for the 
totally irreversible subnetwork of the original "true" chemical reaction net- 
work will, like the SR Graph for the original network, satisfy the requirements 
of Theorem \2.1\ 

This is to say that there is no loss of generality in assuming, for the 
purposes of contradiction, that there is a "true" chemical reaction network, 
containing no reversible reactions, for which the SR Graph satisfies both 
requirements of Theorem 12.11 but for which the fully open extension is dis- 
cordant. 

Hereafter in the proof of Theorem 12.11 then, we assume that all reactions 
in the "true" chemistry are irreversible, that the corresponding SR Graph sat- 
isfies both conditions of Theorem \2.1[ and that, contrary to what the theorem 
asserts, the fully open extension of the true chemistry is discordant. 

Thus there exist fixed a € kerL and nonzero a e R-^ satisfying the re- 
quirements of Definition 14.11 (Recall that the stoichiometric subspace for a 
fully open network is M-^.) R will be understood throughout the proof that all 
references to a and a are relative to this fixed pair, so chosen. 

It will be helpful to divide the proof into subsections: 
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5.1 Preliminaries 

We associate a sign with each species in the following way: A species s ^ 5^ 
is positive, negative, or zero according to whether as is positive, negative, 
or zero. By a signed species we mean one that is either positive or negative. 
Similarly, a reaction y -^ y' ^ M is positive, negative, or zero according to 
whether cxy^yi is positive, negative, or zero. A signed reaction is one that is 
either positive or negative. 

Remark 5.1. Note in particular, that for any "degradation reaction" of the 
kind s -> 0, Definition 14.11 requires that sgn a^^^o = sgn cr^, so sgn s ^ = 
sgn s for every s e 5^ . (A similar situation obtains for any reaction of the 
form ns -^ y, where n is a positive number.) For any "synthesis reaction" 
-^ s. Definition 14. II requires that ao^s = 0, so such reactions are unsigned. 

5.2 The sign-causality graph and causal units 

Hereafter we denote by M* the set of all reactions not of the form s -^ 0, s e S^. 
That is, £$* is the set of reactions that are not degradation reactions. Be- 
cause a is a member of ker L we can write 

E '^y^y'iy' -y)- T, "*-oS = 0. (13) 

Thus, for any particular choice of species s € ^, we have 

E (^y^y'iy's-ys) = (^s^o- (14) 

Now suppose that in (1141) . the species s is positive. From (IT^ and Remark 
15. II we have 

E (^y^y'iys-ys)>0, (15) 

y-*y'eM* 

in which case at least one term on the left must be positive. (At least one 
term on the left is "causal" for the inequalityjj) Terms of this kind can arise 
in precisely two ways: 

(i) There is a positive reaction y -^ y' (i.e., ay^y' is positive) with s e supp?/' 
(so that yg > 0,ys = 0). Recall, however, that for ay^y' to be positive, the 



^Whcn we say that X is causal for the outcome Y, we mean to suggest that X abets 
Y, not necessarily that X, by itself, determines the outcome Y. 
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conditions of Definition 14.11 require that there be a positive species s' in 
suppy. 

In this case, reaction y ^ y' is "causal" for the sign of species s, while 
species s' is "causal" for the sign of reaction y -^ y'. With this in mind, we 
write 



^, y 

S "^ 



y^y' 



y + 



^ s. (16) 



The signs above the species and the reactions indicate their respective signs. 
The complex labels above the "causal" arrows (-^) indicate the complex in 
whose support the adjacent species resides. 

(ii) There is a negative reaction y -^ y' (i.e., ay^y' is negative) with s e supp y 
(so that ys > 0,y'g = 0). But for ay^y' to be negative, the conditions of 
Definition 14.11 require that there be a negative species s' in suppy. As in 
case (i), reaction y ->■ y' is causal for the sign of species s, while species s' is 
causal for the sign of reaction y -> y'. This time we write 



y 



S '^ 



y' 



y + 



"^ 



5. (17) 



On the other hand, suppose that in flT4|) . the species s is negative. From 
f fl4|) and Remark 15.11 we have 

E c^y^y'{ys-ys)<o, (18) 

in which case at least one term on the left must be negative. (At least one 
term on the left is causal for the inequality.) Here again, terms of this kind 
can arise in precisely two ways: 

(i)' There is a negative reaction y -^ y' (i.e., oty^yi is negative) with s e supp y' 
(so that y'g > 0,ys = 0). For ay^y' to be negative, the conditions of Definition 
14.11 require that there be a negative species s' in supp y. Here we write 



/ y 

S '^ 



y' 



y 



(19) 



(ii)' There is a positive reaction y ^ y' (i.e., ay^y' is positive) with s € supp?/ 
(so that ys > 0,y^ = 0). For cty^y' to be positive, the conditions of Definition 
14. II require that there be a positive species s' in suppy. We write 
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+ -<- _ 

s' ^ \y^y'\ ^ s. (20) 



The diagrams drawn in flT6|) . f lT7|l . (IT9|l and (|20|) can be viewed as edge-pairs 
in a bipartite directed graph: 

The sign- causality graph, drawn for the network (relative to the a, a pair 
under consideration) is constructed according to the following prescription: 
The vertices are the signed species and signed (non-degradation) reactions. 
An edge -^ is drawn from a signed species s' to a signed reaction y -^ y' 
whenever s' is contained in supp y and the two signs agree; the edge is then 
labeled with the complex y. An edge -^ is drawn from a signed reaction 
y -^ y' to a signed species s in either of the following situations: (i) s is 
contained in supp y' and the sign of s agrees with the sign of the reaction; in 
this case the edge carries the label y' or (ii) s is contained in supp y and the 
sign of s disagrees with the sign of the reaction; in this case the edge carries 
the label y. It is understood that the signed species and the signed reactions 
are labeled by their corresponding signs. 

By a causal unit we mean a directed two-edge subgraph of the sign- 
causality graph of the kind s' ^ R^ s, where R denotes a reaction. (We will 
often designate a reaction by the symbol R when there is no need to indicate 
the reactant and product complexes.) The species s' is the initiator of the 
causal unit s' -^ R -^ s, while s is its terminator. It is not difficult to see that 
the initiator and terminator of a causal unit must be distinct species|l| 

Causal units are of the four varieties shown in f lTB]) . f lT7|) . flTI?]) and fl2U]) . 
Of these, (TT6|) and (fT9|) carry distinct complex labels on the two edges. On 
the other hand, (TT7|) and ( 120|) carry identical labels on the edges. By a c-pair 
causal unit we mean a causal unit of the second kind. (As with the Species- 
Reaction Graph, the term is meant to be mnemonic for complex pair.) 

Remark 5.2. It is important to note that a c-pair causal unit results in a 
change of sign as the edges are traversed from the initiator species to the 
terminator species. Otherwise, a causal unit is characterized by retention of 
the sign. 

Remark 5.3. It should be clear that, apart from the direction -* imparted to 
its edges, the sign-causality graph corresponding to a and a can be identified 



^Recall that, by supposition, a species can appear on only one side of a reaction. 
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with a subgraph of the undirected Species-Reaction Graph (which we will 
sometimes refer to as the sign-causality graph's "counterpart" in the Species- 
Reaction Graph). Moreover, every fixed-direction edge-pair in the SR Graph 
has direction consistent with that imparted by the -^-relation (because each 
proceeds from a reactant species to an irreversible reaction to a product 
species). 

5.3 The sign-causality graph must contain a directed 
cycle, and all of its directed cycles are even. 

By supposition a is not zero, so there is at least one signed species, say Si. 
From the discussion in Section 15.21 it is clear that Si is the terminator of a 
causal unit S2 ^ R2 ^ si, where S2 is distinct from the initiator si. Because 
S2 is also signed, it too must be the terminator of a causal unit S3 -* i?3 -^ S2, 
where the signed species S3 is distinct from S2. Continuing in this way, we 
can see that there is a directed sequence of the form 

■■■ -^ Sn-^ Rn '^ Sn-1 "^ Rn-V -^ S3 -> i?3 "^ S2 "^ i?2 '^ Si, (21) 

with Si 4^ Sj+i. 

Because the number of species is finite, two non-consecutive species in the 
sequence must in fact coincide, which is to say that the sign- causality graph 
must contain a directed cycle. Moreover, it is easy to see that each vertex in 
the sign- causality graph resides in a cycle or else there is a cycle -^-upstream 
from it. 

With this in mind, we record here some vocabulary that will be useful 
in the next section: A source is a strong component of the sign-causality 
graph whose vertices have no incoming edges originating at vertices outside 
that strong component. Clearly, every component of the sign-causality graph 
contains a source, and every vertex in a source resides in a directed cycle. 

As with the SR Graph, we say that a (not necessarily directed) cycle in 
the sign-causality graph is even if it contains an even number of c-pairs. 

Lemma 5.4. A (not necessarily directed) cycle in the sign- causality graph 
that is the union of causal units is even. In particular, every directed cycle 
in the sign- causality graph is even. 

Proof. If we traverse the cycle beginning at a species s° and count the number 
of species-sign changes when we have returned to s°, that number clearly 
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must be even. But, if the cycle is the union of causal units, the number 
of sign changes is identical to the number of c-pair causal units the cycle 
contains (Remark I5.2p . Clearly, every directed cycle in the sign-causality 
graph is the union of causal units. D 

Because a directed cycle in the sign-causality graph is even, its (orient- 
table) cycle counterpart in the undirected Species-Reaction Graph (Remark 
15. 3p will also have an even number of c-pairs. Because the sign-causality 
graph must contain a source, which in turn must contain a directed cycle, 
we now know that a reaction network is concordant if its Species- Reaction 
Graph contains no orientable even cycles. 

Remark 5.5. In fact, a source in the sign- causality graph, when viewed in the 
SR Graph, must be a critical subgraph. That this is so follows from the fact 
that a source is a strong component of the sign-causality graph and therefore 
is the union of -^-directed cycles. These even cycles, viewed in the SR Graph, 
have a consistent orientation, the orientation conferred by the directed-cycle 
"^-orientation in the sign-causality graph, which in turn is consistent with 
any fixed-direction edge-pairs the SR Graph might contain. 

In the next section we begin to consider what happens when the Species- 
Reaction Graph does contain at least one orientable even cycle (and therefore 
at least one critical subgraph). We will want to show that if the fully open 
network under consideration satisfies the hypothesis of Theorem 12.11 then 
the very existence of a source in the putative sign-causality graph becomes 
impossible Ifl 



5.4 Inequalities associated with a source; stoichiomet- 
ric coefficients associated with edges in the sign- 
causaUty graph 

Consider a source in the sign-causality graph having J^q c ^ as its set of 
species nodes and ^o "^ ^* as its set of reaction nodes. If s is a positive 
species in the source, then 015p can be written 



*As a matter of vocabulary, a strong component of a directed graph — in particular a 
source — is, formally, a set of vertices in the graph, but here, less formally, we will also 
associate the source with the subgraph of the sign-causality graph obtained by joining the 
vertices of the source with directed edges inherited from the parent graph. 
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E Oiy^y'iy's-ys) + E O!y^y'iys-ys)>0. (22) 

Now suppose that a term in the second sum on the left, corresponding to 
reaction y -^ y', is not zero. Because y -^ y' is not a vertex of the source, any 
edge of the sign-causahty graph that connects species s to reaction y ^ y' 
must point away from s. Thus, the reaction y -^ y' cannot be causal for the 
positive sign of species s, so the term 

(^y~*y'{y's-ys) 
must be negative. This implies that (12 2 p can obtain only if we have 

^ ay^y,{y',-y,)>Q. (23) 

When s is a negative species in the source, we can reason similarly to write 

E ay^y'Ws-ys)<Q- (24) 



Remark 5.6. Note that, for a particular s e S^q a nonzero term in (123|) or (12] 
might not correspond to any edge of the sign-causality graph (as when, for a 
particular reaction y ^ ?/', s is a member of supp y' and disagrees in sign with 
ay^y>). If s is a positive species, then the term in question must be negative, 
and hence it can be removed from (123|) without changing the sense of that 
inequality. Similarly, if s is a negative species, the term in question is positive 
and can be removed from fl2^ without changing the inequality's sense. In 
what follows below, we shall assume that such removals have been made, so 
that every term in (125]) or ([21]) corresponds to an edge in the sign- causality 
graph. 

Recall that a directed edge of the sign-causality graph from a species s 
to a reaction y ^ y' is always of the form 



y 



y^y' ■ 



That is, the edge-label of such a species-to-reaction edge is invariably the 
reactant complex, y. Note that species s has a positive stoichiometric coef- 
ficient, ys, in that complex. On the other hand, reaction-to-species edges of 
the sign- causality graph are of two kinds: 

1 y I 1 y' 

y ^ y' ^ s or y -^ y' ^ s. 
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In either case, the species s has a positive stoichiometric coefficient (either 
Us or y'g) in the edge-labehng complex. 

Hereafter, for a species-to-reaction edge s -^ R of the sign- causality graph 
we denote by Cs^r the (positive) stoichiometric coefficient of species s in the 
corresponding edge-labeling complex. For a reaction-to- species edge R -^ s 
we denote by fji^s the (positive) stoichiometric coefficient of species s in its 
edge-labeling complex. 

For the sign-causahty graph source under consideration, we will in fact 
need a small amount of additional notation: For each species s € ^q "we 
denote by ^o ^ s the set of all edges of the source that are incoming to 
s and by s -> ^q the set of all edges of the source that are outgoing from 
s. In light of Remark 15.61 and in view of notation we now have available, 
some analysis will indicate that the inequalities given by fl23|) and f l24|) can 
be written as a single system: 



X! fn-slaRl - X! ^s-^nlanl > 0, Vs e ^o- 



(25) 



Our aim will be to show that when the conditions of Theorem 12. II obtain, 
the inequality system (125|) cannot be satisfied. The following remark will 
play an important role. 

Remark 5.7. Consider a source in the sign-causality graph, and suppose that 

Si -> -Ri -^ S2 -^ -R2 ■ ■ ■ '^ ■§„ -^ i?„ -> Si (26) 

is one of its -^-directed (and therefore even) cycles. Note that the -^-orientation 
of the cycle also provides an orientation of the cycle's counterpart in the 
Species-Reaction Graph, for it is consistent with directions of the fixed- 
direction edge-pairs in the SR Graph. (Remark 15. 3p 

Thus, when condition (i) of Theorem \2.1\ is satisfied, the stoichiometric 
coefficients of the directed cycle ( l26|) must satisfy the condition 

^si--^Ri^S2-^R2 • ■ • ^s„-^R„ 
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5.5 The decomposition of a source into its blocks; the 
block-tree of a source 

Here we draw on graph theory vocabulary that is more-or-less standard ^ jj 
A separation of a connected graph is a decomposition of the graph into two 
edge-disjoint connected subgraphs, each having at least one edge, such that 
the two subgraphs have just one vertex in common. If a connected graph 
admits a separation (in which case it is separable), then the common vertex 
of the separation is called a separating vertex of the graph. A graph is 
nonseparable if it is connected and has no separating vertices. A maximal 
nonseparable subgraph of a graph is called a block of the graph. In rough 
terms, a connected graph is made up of its blocks, pinned together at the 
graph's separating vertices. 




Figure 2: A hypothetical source with five blocks 



^Our focus is on bipartite graphs, so in the discussion of terminology it will be under- 
stood that no vertex is adjacent to itself via a self- loop edge. 
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A source in the sign-causality graph can be decomposed into its blocks, 
which we call source-blocks. Because a source is strongly connected, each of its 
blocks is strongly connected (and non-separable). To illustrate some of these 
ideas we show in Figure |2] a hypothetical source with five blocks. (Although 
they are inconsequential to the block decomposition, the arrows in the figure 
are meant to connote the -^-relation.) The separating vertices in the figure 
are R2, S3, S4, and 5*7. An example of a block is the subgraph having vertices 
{Si,S2,Ri,R2,R3} and edges {Si -> i?2,-R2 -^ S2,S2 -^ Ri,Ri ^ Si, Si ^ 

R31R3 '^ 'S'2}- 

The block-tree of a connected graph depicts the way in which the various 
blocks of the graph are joined at its separating vertices. More precisely, the 
block-tree |1] of a connected graph is a bipartite graph whose vertices are of 
two kinds: symbols for the graph's blocks and symbols for the graph's sep- 
arating vertices. If a block contains a particular separating vertex, then an 
edge is drawn from the block's symbol to the symbol for that separating ver- 
tex. In Figure [3] we show the block-tree for the hypothetical source depicted 
in Figure [2J For reasons that will be made clear later on, we have denoted 
the five blocks in the block-tree by the symbols RB1,SB1,SB2,RB2, and 
RB3. (We note that Figures [2] and [3] are somewhat unrepresentative, for 
it might happen that three or more source-blocks are adjacent to the same 
separating vertex.) 




Figure 3: The block-tree for Figure [2] 

A block is an end block if it contains no more than one separating vertex 
of the original graph. In this case, the block's symbol is a leaf of the block- 
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tree graph. The end-blocks (leaves) in our example are those represented by 
RB1,SB2, and RB3. 

5.6 Properties of source-blocks 

Our aim in this subsection is to show that when condition (ii) of Theorem 12 .11 
is satisfied, the internal structure of source-blocks must have a certain degree 
of simplicity. Condition (ii) will exert itself primarily through the following 
proposition, which is proved in Appendix Rl 



Proposition 5.8. Suppose that, for the reaction network under consider- 
ation, the Species-Reaction Graph satisfies condition (ii) of Theorem \2.1[ 
Then, in any source-block of the sign- causality graph, at most one of the 
following can obtain: 

(i) There is a reaction vertex having more than two adjacent species ver- 
tices. 

(ii) There is a species vertex having more than two adjacent reaction 
vertices. 

The proposition provides motivation for the following definition: 

Definition 5.9. A source-block in the sign-causality graph is a species block 
(S-block) if each species node is adjacent to precisely two reaction nodes. A 
source-block in the sign- causality graph is a reaction block (R-block) if each 
reaction node is adjacent to precisely two species nodes. 

Proposition l5.8l tells us that when condition (ii) of Theorem 12. 6 1 is satisfied, 
every block within the sign-causality graph source is either an S-block or an 
R-block (or both in the case that the source-block is simply a single cycle). 
Note that in Figure [3] we have labeled the hypothetical source-blocks as RBn 
or SBn according to whether the corresponding source-block in Figure [2] is 
an R-block or an S-block. 

We conclude this section with two more propositions. Neither is essential 
to the proofs of the main theorems of this paper, but they provide some 
additional and not-so-obvious properties of a sign- causality graph source. 

The following proposition, proved in Appendix \^ does not presuppose 
that condition (ii) of Theorem 12.11 is satisfied. Rather, it tells us about 
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properties of R-blocks or S-blocks that might exist within a sign-causahty- 
graph source. We aheady know that every directed cycle is even, but the 
proposition tells us that all cycles within a source's R-blocks and S-blocks 
are even. 

Proposition 5.10. Every (not necessarily directed) cycle that lies within an 
S-block or an R-block of a sign- causality graph source is even. 

The following proposition is a direct consequence of the two preceding 
ones: 

Proposition 5.11. Suppose that, for the reaction network under consider- 
ation, the Species- Reaction Graph satisfies condition (ii) of Theorem \2.1\ 
Then in any source of the sign- causality graph every cycle is even. 

5.7 Properties of an end S-block 

Recall that an end S-block in a source is an S-block that contains at most 
one separating vertex of the source. We consider properties of such an end S- 
block, designated ESB. There are three mutually exclusive possibilities: ESB 
contains no separating vertex at all; it contains just one separating vertex, 
and it is a reaction vertex; or it contains just one separating vertex, and it is 
a species vertex. For our purposes the first two possibilities can be treated 
together, while the third requires other considerations. 

In fact, we show that when condition (i) of Theorem 12.11 holds, the first 
two possibilities cannot obtain; moreover, if the third obtains, we get sharp- 
ened information about the inequality in ( 12^ corresponding to species at the 
block's separating vertex. 

5.7.1 Possibilities 1 and 2: ESB contains no separating vertex of 
the source or it contains a separating reaction vertex of the 
source 

Because it is strongly connected, ESB must contain a directed (and even) 
cycle, which we take to be 

Si ^ Rl ^ $2 ^ R2- ■ ■ ^ Sn ^ Rn ^ Si- (28) 

Because ESB is a species-block and because the block contains no separating 
species-vertex of the source, each of the species in the block (and in the 
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chosen cycle) is adjacent to at most two reactions of the block. Thus, the 
inequalities f l^ corresponding to Si, S2, . • . , s„ reduce to 

/i?„-.si|tt/?„| -Csi^Rja^J > 0, 
fR^^S2\aR^\-es2-^R^\aR^\ > 0, (29) 

(Recall that, for a species-to- reaction edge s --^ i? of the sign- causality graph 
we denote by Cs^r the stoichiometric coefficient of species s in the corre- 
sponding edge-labeling complex. For a reaction-to-species edge /? -> s we 
denote by //j^s the stoichiometric coefficient of species s in its edge-labeling 
complex.) By sequentially invoking these inequalities from top to bottom, 
we can deduce from the last of them that 



fRl'^S2fR2~'R3 ■ ■ ■ fR-n- 



■si 



^Sl~rR\(^S2'^R2 ■ ■ ■ ^Sn-^Rn 



1 |aijj>0. (30) 



However, when condition (i) of Theorem 12.11 holds, ( l30l) cannot obtain: By 
supposition |q;r„| is positive. Given the -^-orientation in the SR Graph, 
the even cycle ( 128|) (viewed in the SR Graph) cannot be stoichiometrically 
expansive, so the first factor on the left of (!30|) is either zero or negative. 
(Recall (127|) .) Thus, we have a contradiction of (!30|) . 

We conclude, then, that when condition (i) of Theorem 12. II holds, an end 
S-block must contain a separating species vertex of the source. We investigate 
next what happens in that case. 

5.7.2 Possibility 3: ESB contains a separating species vertex of 
the source 

Here again we consider a directed cycle in ESB, labeled as in (l28ll . If ESB's 
(unique) separating species vertex is not in the cycle, then we would again 
obtain a contradiction, just as in § 15.7.11 We suppose, then, that species s„ 
is the separating vertex. Because ESB is a species-block, all other species 
vertices of the cycle are adjacent to precisely two reaction vertices, and those 
are in the cycle. Thus, all inequalities but the last in f l29|) remain unchanged. 
On the other hand, s„ is adjacent not only to Rn-i and i?„ but also to 
certain reactions from nearby blocks sharing s„ as a species. Recall that Mq 
is the set of all reactions in the source under study. We denote by S^esb the 
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set of reactions in ESB and by ^q \ ^esb the set of all source reactions not 
in ESB. Moreover, we let Mq \^esb ^ ■§„ and s„ ^ Mq \Mesb be the sets of 
source edges residing outside of ESB that are, respectively, directed toward 
and away from species s„. 

In this case, the last inequality in (|29|) must be replaced by 



(31) 
Instead of (!30!) . this time we deduce the inequality 

(J R\~^S2J R2-~^Ri • ■ • JRn^Sl 1 I |„ I , /'QO^ 
1 FRnI + (32) 

When condition (i) of Theorem [2]T] obtains, however, the first term on the left 
cannot be positive for reasons given in §5.7.11 Thus, we arrive at the following 
inequality, which relates entirely to source edges disjoint from ESB: 

E fR~'s„\oiR\ - Y. es„^ij|aK| > 0. (33) 

For species Sn this amounts to a sharpened form of its counterpart in ( l25l) . a 
form we will draw upon later on. 

5.8 Properties of an end R-block 

We begin this subsection with an important proposition about R-blocks. A 
proof is provided in Appendix [Bl 

Proposition 5.12. Suppose that, in a sign- causality graph source, an R- 
block has species set 5^* , and suppose that no directed cycle in the block is 
stoichiometrically expansive relative to the ^ orientation. Then there is a 
set of positive numbers {Ms}s^y'* such that, for each causal unit s -^ R-^ s' 
in the block, the following relation is satisfied: 

fR^s'Ms' - Cs^rMs < 0. (34) 
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Now we consider an end R-block in the source under consideration, des- 
ignated ERB. We denote by 5^* the set of species in ERB. In consideration 
of the source inequahty system fl25l) . we restrict our attention to just those 
inequahties corresponding to species in 5^*: 

E /i?-s|aij| - X! e,^i?|ai?l > 0, Vs e ^*. (35) 

We suppose that condition (i) of Theorem 12.11 is satisfied, so, by virtue of 
Remark [521 ^^ can choose {Ms}sej^* as in Proposition 15.121 If we muhiply 
each inequahty in ([25]) by the corresponding M^ and sum, we get the single 
inequality shown in (I5B1) . 



e(e/ 



,j^,M,|a^| - E M,e,^R\aR\ | > (36) 



As in §5.71 there are three possibilities: ERB contains no separating ver- 
tex; it contains just one separating vertex, and it is a reaction vertex; or 
it contains just one separating vertex, and it is a species vertex. We will 
show that neither of the first two possibilities can obtain. Then, as in §5.7.2[ 
we will show that, if the third possibility is realized, the inequality in ( 125|) 
corresponding to the species at the separating vertex can be sharpened to 
considerable advantage. 

5.8.1 Possibilities 1 and 2: ERB contains no separating vertex of 
the source or it contains a separating reaction vertex of the 
source 

Because ERB is an R-block, each reaction is adjacent to precisely two species 
in the set ^*, which is to say that each reaction in ERB sits at the center of 
precisely one causal unit in ERB. (Recall that ERB is strongly connected.) 
Moreover, if either of the first two possibilities should obtain, no species is 
adjacent to a reaction not in ERB. In these cases, the inequality (I36p can be 
rewritten. Let '^ be the set of causal units within ERB. Then (I36|) can be 
made to take the form shown in ([3j 



E (/ii^s'M,,-e,^HM,)M > (37) 

However, ([37]) is contradicted by the attributes given to the set {M^g^y* in 
Proposition I5.12[ 
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Thus, neither of the first two possibihties can obtain. ERB must contain 
a separating species vertex of the source. We examine next what can be said 
in that case. 

5.8.2 Possibility 3: ERB contains a separating species vertex of 
the source 

Suppose that ERB contains a species vertex s* that is a separating vertex of 
the source under consideration. When condition (i) of Theorem 12.11 obtains, 
we can again choose positive numbers {Ms}se,y* to satisfy the requirements 
of Proposition 15. 12^ and the inequahty fl5^ remains in force. On the other 
hand the passage from fl5Bl) to fl57|) becomes confounded by the fact that 
species vertex s* is now adjacent to source edges not residing in ERB. We 
denote by Mere the set of reactions in ERB and by Mq \ Mere the set 
of all source reactions not in ERB. Moreover, we let .^o ^ ^erb ^ s* and 
s* ^ ^0 ^^ERE be the sets of source edges residing outside of ERB that are, 
respectively, directed toward and away from species s*. As before, we let "^ 
be the set of causal units in ERB. In this case, f p6|) can be recast as f p8|l . 

Y, ifR-^s'Ms' - es-^RMs)\aR\ + 

(38) 

Ms*{ Yj fR-^s*\0CR\ - Yj (is*-^R\c^R\) > 

Recall that the set {Ms^s^y* was chosen to satisfy the requirements of 
Proposition I5.12[ so the first sum in f p8|) cannot be positive. Then, because 
Ms* is positive, we must have 

Y fR-sAoinl - Y (is*'^R\oiR\ > 0. (39) 

Note that this is a strengthened counterpart of the inequality in (12 5 p corre- 
sponding to species s*, a counterpart that makes reference only to reactions 
external to the end reaction block ERB. 



5.9 The concluding argument: Leaf removal 

We begin this subsection with a review of what was established in §5.71 and 
§5.81 An end block in a sign-causality graph source, be it an S-block or an 
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R-block, must contain a species separating vertex of the source. This imphes 
that the hypothetical source depicted in Figure [2] (and having block-tree 
depicted in Figure [3]) cannot in fact be a source, for it has an end block, 
corresponding to RBI in Figure [3l that does not contain a separating species 
vertex of the putative source. (Each of the remaining end blocks does contain 
a separating species vertex.) 

Moreover, we have established that the inequality in 025 p corresponding 
to a separating species vertex s* in an end block EB, be it an S-block or an 
R-block, can be strengthened to a form that makes no mention of reaction 
inside EB: 

Yj /i?-s*|aij| - Yj (is*-,R\oiR\ > 0. (40) 

Here ^eb is the set of reactions in EB, and ^o \ Meb the set of source 
reactions not in EB. Moreover, ^o ^ ^eb ^ s* and s* -* 3?o \ ^eb are the 
sets of source edges residing outside of EB that are, respectively, directed 
toward and away from species s*. 

Now if EB is an end block in the source under consideration, we can 
replace the inequality in (!25|) corresponding to the unique separating species 
vertex s* in EB with its strengthened form shown in (HOl) . Thereafter, we 
can restrict the now-modified inequality system (1251) to just those inequalities 
corresponding to s* and to species residing outside EB. In effect, the resulting 
reduced system of inequalities corresponds to a subgraph of the source with 
the end block EB removed, but with s* retained. Viewed in the source's block 
tree, this corresponds to removal of a leaf along with that leaf's adjacent 
species separating vertex (when that separating vertex is not adjacent to a 
different leaf). 

It is not difficult to see that the arguments in §5.71 or §5.81 can then be 
applied to any end block EB' of the resulting source subgraph to produce 
a still smaller but strengthened inequality system corresponding to a still 
smaller source subgraph, a subgraph resulting from removal of EB'. 

The process can be continued, amounting to a sequential pruning of the 
source's block-tree, with each stage corresponding to removal of a (perhaps 
new) leaf and perhaps its adjacent species separating vertex. The process 
will terminate when only one block remains, with the final source subgraph 
having no separating vertex at all. In this case, by arguments of §5.7.11 or 
§5.8.11 the corresponding reduced inequality system cannot be satisfied, and 
we have a contradiction. 
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This completes the proof of Theorem 12. 1[ 



6 Proof of Theorem 2.6 



Here we prove Theorem I2.6[ which is repeated below: 

Theorem 12.61 A fully open reaction network is strongly concordant if its 
true- chemistry Species- Reaction Graph has the following properties: 

(i) Every even cycle is an s-cycle. 

(a) No two even cycles have a species-to-reaction intersection. 

We begin by recalling the definition of strong concordance for an arbitrary 
reaction network {^,'ia',M} , not necessarily fully open, with stoichiometric 
subspace S c M.-^ . Again we let L : M-^ ^ S" be the linear map defined by 

La= Y. Oiy^y'(y'-y)- (41) 

Definition 6.1. A reaction network {=5^, ^,^} with stoichiometric subspace 
S is strongly concordant if there do not exist a e ker L and a non-zero a e S 
having the following properties: 

(i) For each y -^ y' such that ay^y' > there exists a species s e supp {y-y') 
for which sgn cr^ = sgn (y - y')s. 

(ii) For each y ^ y' such that ay^y' < there exists a species s e supp {y-y') 
for which sgno"s = -sgn (y - y')s. 

(iii) For each y -> y' such that ay-^y' = 0, either (a) as = for all s e supp y, or 
(b) there exist species s, s' e supp {y - y') for which sgn ag - sgn {y - y')s 
and sgncT^. = -sgn {y - y')^'. 

The following lemma will take us just short of a proof of Theorem 12.61 
Proof of the lemma is given in Appendix O 

Lemma 6.2. Suppose that a fully open reaction network with true chemistry 
{S^j'if,^} is not strongly concordant. Then there is another true chemistry 
{^,'t^,M} whose fully open extension is discordant and whose SR Graph is 
identical to a subgraph of the SR Graph for {^,'^,^} , apart perhaps from 
changes in certain arrow directions within the reaction vertices. 

41 



Proof of Theorem 12.61 proceeds from Lemma 16.21 in the following way: 

Proof of Theorem \2.(A Suppose that the SR Graph for a true chemistry 
{,y,^,^} satisfies conditions (i) and (ii) of the theorem but that, contrary 
to what is to be proved, the fully open extension of {^, ^, ^} is not strongly 
concordant. Note that when the SR Graph of a true chemistry {^,^,=^} 
satisfies conditions (i) and (ii) of Theorem 12. 6[ so will any subgraph of that 
SR Graph. On the other hand, neither of those conditions depends upon the 
direction of arrows in the reaction nodes. Thus, the SR Graph of the true 
chemistry l^,*^,^} given by the lemma will also satisfy conditions (i) and 
(ii). But then, by Corollary 12.41 of Theorem l2.lt the fully open extension of 
{y,^^M} could not be discordant, and we have a contradiction. D 

7 Extensions of the main theorems to net- 
works that are not fully open 

It is the purpose of this Section to elaborate on remarks made in §1.2.21 

When the SR Graph drawn for a true chemical reaction network satisfies 
the hypotheses of Theorems 12.11 or 12. 6[ these theorems tell us that the net- 
work's fully open extension is concordant (or strongly concordant). In this 
case, the fully open network inherits the many attributes described in [18] 
that accrue to all concordant (or strongly concordant) networks. We would 
like to know circumstances under which these theorems can be extended in 
their range to give concordance information about networks that are not fully 
open. 

More generally, we would like to know conditions under which, for a given 
network, concordance or strong concordance of the network's fully open ex- 
tension implies concordance or strong concordance of the network itself]^ 
This is a question quite separate from SR Graph considerations. However, 
when the network satisfies such conditions and its underlying true chemistry 
SR Graph satisfies the hypothesis of Theorems 12.11 or 12. 6[ then the concor- 
dance properties ensured by those theorems for the fully open extension will 
be inherited by the original network. 

In [18] we showed that a normal network is concordant (strongly concor- 
dant) if its fully open extension is concordant (strongly concordant). Nor- 



""^"We do not preclude the possibility that the original network contain some degradation 
reactions of the form s -*■ 0. 
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mality is a mild structural condition given in Definition 17.11 below. In [S] it 
was shown that every weakly reversible network is normal. 

Theorems 17.101 and 17.111 below tell us that these same results also obtain 
for the still larger class of weakly normal networks. (See Definition 17.31 ) In 
particular, a weakly normal network is concordant (strongly concordant) if 
its fully open extension is concordant (strongly concordant). 

This improvement on results in [TS] is helpful in itself, but it also has 
significance in another direction: We show in §7.11 that the class of weakly 
normal networks is synonymous with the very broad class of nondegenerate 
networks (Definition 17. 6p . which was described in the Introduction. Thus, 
any nondegenerate reaction network with a concordant (strongly concordant) 
fully open extension is itself concordant (strongly concordant). Moreover, we 
also show that, with respect to the possibility of concordance, degenerate 
reaction networks are not worth considering, for they are never concordant. 

In §7.31 we provide computational tests that serve to affirm network nor- 
mality and weak normality (or, equivalently nondegeneracy) . 

7.1 Network normality, weak normality, and nonde- 
generacy 

Definition 7.1. Consider a reaction network {^,'^,^} with stoichiometric 
subspace S. The network is normal if there are q € Mf^ and i] e M-f such that 
the linear transformation T : S -^ S defined by 

T(r:= E VyMy*^)iy'-y) (42) 

is nonsingular, where "*" is the scalar product in M-^ defined by 

x*x' := Y, ^sa^sa^s- (43) 

Remark 7.2. As indicated earlier, it was shown in P] that every weakly 
reversible network is normal. Reference j8j also contains structural conditions 
that ensure normality for certain "partially open" networks that are not 
weakly reversible. 

In preparation for the next definition we note that (H2]) can be written as 

Ta:= E ivy^y'y°<i)-^iy'-y), (44) 
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where "•" indicates the standard scalar product in 

Definition 7.3. Consider a reaction network {^5^,^,^} with stoichiometric 
subspace S. The network is weakly normal if, for each reaction y -^ y', 

there is a vector Py^y' e M+ with snpp Py^yr = suppy such that the hnear 
transformation T : S ^ S defined by 

T(t:- Y, Py^y<y'-y) (45) 

is nonsingular. Here "•" is the standard scalar product in M'^. 

Remark 7.4. A reaction network that is normal is also weakly normal. In 
fact, if g e Mf and r] e Mf satisfy the requirements of Definition I7.H then 
the choice Py^y' '■= Vy^y'V ° q,^y -^ y' ^ ^ will satisfy the requirements of 
Definition 17.31 On the other hand, a weakly normal network need not be 
normal. An example is given by 

C^ A + B^D^2A, (46) 

which is weakly normal but not normal. Because every weakly reversible 
network is normal, it follows that every weakly reversible network is also 
weakly normal. 

Remark 7.5. For readers familiar with standard language of chemical reaction 
network theory p^, a network cannot be normal if 

t-£-6 > 0, (47) 

where t is the number of terminal strong linkage classes, i is the number of 
linkage classes, and 6 is the deficiency. This follows without much difficulty 
from [12]; see also [10] • For network ( 146|) . t = 2, i = 1, and 6 = 0, so normality 
is precluded by condition ( H71) . Network ( H6l) illustrates, however, that the 
same condition does not also preclude weak normality. 

Definition 7.6. A reaction network is nondegenerate if there exists for it 
different iably monotonic kinetics ( §3.2p such that at some positive composi- 
tion c* the derivative of the species-formation-rate function df{c*) '■ S -^ S is 
nonsingular. Otherwise, the network is degenerate. 
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Proposition 7.7. A reaction network is nondegenerate if and only if it is 
weakly normal. 

Proof. Suppose that a network {S^,^,^} is nondegenerate. Then there 
is for the network a kinetics J^ such that, at some composition c* € Mf^, 
the kinetics is differentiably monotonic and, moreover, the derivative of the 
species-formation-rate function, df(c*) ■■ S -^ S, is nonsingular. In this case, 
for each a € S 

df{c*)a= Y. VJ(ry^y>{c*)-a{y'-y), (48) 

where the components of V^^y'(c*) have the non-negativity properties that 
follow from different iable monotonicity ( §3.2p . By taking 

Py^y' = V^Xj^y' (c* ) , ^y^y'^M 

we can establish that the network is weakly normal. 

On the other hand, suppose that the network {^, ^, .^j is weakly normal 
and, in particular, that the set {Py^y'}y^y'(iy2 satisfies the requirements of 
Definition 17.31 Let J(f be the (differentiably monotonic) kinetics defined by 

jey^y,{c) := (fy-y' , ^y^y'^M, 

and let c* e Mf^ be such that c* = 1 for each s € ^. Note that 

VJfy^y'(c* ) = py^y, , ^y -^y' e^. 

From this, fHHj) . and the properties of the set {Py^y'}y^y'eM given by Definition 
17.31 it follows that df{c*) '■ S -^ S, is nonsingular, whereupon the network is 
nondegenerate. D 

Remark 7.8. Note that in the proof that nondegeneracy implies weak nor- 
mality we did not actually require that the network be nondegenerate. In 
particular, we did not require that the kinetics J^ be differentiably mono- 
tonic at all positive compositions, only that it be differentiably monotonic at 
one composition, c* (and, of course, that df{c*) be nonsingular). However, 
when these apparently milder conditions are invoked, the network must be 
nondegenerate nevertheless: The seemingly milder conditions result in weak 
normality, and, as the second part of the proof indicates, weak normality 
implies nondegeneracy. 
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The following proposition indicates that a network that is not weakly 
normal (or, equivalently, is degenerate) has no chance of being concordant. 

Proposition 7.9. A reaction network that is not weakly normal is discor- 
dant. Equivalently, every degenerate network is discordant. 

Proof. Suppose that a reaction network {^,'t^,M} is not weakly normal 
(and, in particular, is not normal). From Definition 17.31 it follows that, for 
the special choice Py^y' = y, ^y -^ y' e ^, the corresponding map T : S ^ S 
given by (H5|) must be singular. This is to say that there is a nonzero a* e S 
such that 

E ycr*{y'-y) = 0. 

Now let a e M.^ be defined by cty^y' = y ■ a*, ^y ^ y' e ^. Then, in view of 
Definition 14. H the pair consisting of a and a* serve to establish discordance 
of the network under consideration. D 

7.2 Concordance of a network and of its fully open 
extension 

The following theorems about weakly normal networks amount to straight- 
forward extensions of theorems in [18] about normal networks; the proofs are 
almost identical. Although these theorems give the former ones a somewhat 
greater range, their main interest lies in the fact that they can be stated in 
terms of the more tangible, but equivalent, notion of network nondegeneracy. 

Theorem 7.10. A weakly normal (or, equivalently, nondegenerate) network 
is concordant if its fully open extension is concordant. In particular, a weakly 
reversible network is concordant if its fully open extension is concordant. 

Theorem 7.11. A weakly normal (or, equivalently, nondegenerate) network 
is strongly concordant if its fully open extension is strongly concordant. In 
particular, a weakly reversible network is strongly concordant if its fully open 
extension is strongly concordant. 

Taken together. Proposition 17. 7[ Proposition 17. 9[ Theorem 17. 10[ and The- 
orem 17.111 tell us that, in the class of networks that have a concordant 
(strongly concordant) fully open extension, the concordant (strongly con- 
cordant) ones are precisely the nondegenerate ones: 

46 



Corollary 7.12. Consider a reaction network that has a concordant (strongly 
concordant) fully open extension. Then the original network is concordant 
(strongly concordant) if and only if it is nondegenerate. 

Remark 7.13. It is a consequence of Theorem 17.101 and results in [18] that the 
dynamical statements (i) - (iii) of Theorem 12.81 hold true for any nondegen- 
erate reaction network whose fully open extension is concordant, not merely 
those that satisfy the SR Graph conditions of Theorem 12.11 In particular, 
the "all or nothing" observation of Remark 12.91 still obtains. The SR Graph 
conditions of Theorem 12.11 merely suffice to ensure concordance of the fully 
open extension. 

7.3 Tests for network normality, weak normality, and 
nondegeneracy 

Here we provide some computational means to affirm normality and weak 
normality (or, equivalently, nondegeneracy) of a reaction network. Recall 
that the rank of a network is the rank of its set of reaction vectors. 

Proposition 7.14. A reaction network {^,^,,^} of rank r is weakly normal 

(or, equivalently, nondegenerate) if there is a set of r reactions {yi -^ yl}i=i...r 

y 

and a set of vectors {pi}i=i...r <= '^+ with snpppi = suppy,, i = l...r, such 

that the matrix 

h•(^^%■)],,.l.... (49) 

has nonzero determinant. 

Proof of the proposition is provided at the end of this subsection. For 
a particular choice of r reactions one can readily construct the matrix 
in terms of symbols for the species-wise components of the vectors {pi}i=i 
and one can then calculate the determinant of the matrix as a polynomial 
in those same symbols. If, for even one choice of r reactions, the resulting 
determinant is not identically zero, then the network is weakly normal. 

We remark in passing that a nonzero determinant will require that the 
r reaction vectors {yl - yi}i=i...r be linearly independent. Of course, such 
independent reaction vector sets will invariably exist for a network of rank r. 

A special choice of {pi}i=i...r, one that obviates the need for symbolic 
computation, is invoked in the following corollary. This choice will often 
suffice to establish weak normality. In fact, when the condition in Corollary 
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17.151 below is satisfied, the network will not only be weakly normal but also 
normal. (See Remark 17. 181 following the proof of Proposition 17. 141 ) 

Corollary 7.15. A reaction network {^,^,=^} of rank r is weakly normal 
(and, in fact, normal) if there is a set of r reactions {yi -^ y'j}i=i...r such that 
the matrix 

h-^'rvA,.,... (50) 

has nonzero determinant. 

Example 7.16. Here we apply Corollary 17.151 to network ([T]), the rank of 
which is 4. For the four-reaction set 

{P^A + B,Q^B + C,2A^C,C + D^Q + E}, (51) 

it is easy to calculate that the matrix in the corollary has a determinant of 
4, so the network is weakly normal. Recall that Theorem 12 . 1 1 established the 
concordance of the fully open extension of network ([T]) . Because network ([T]) 
is weakly normal, concordance of its fully open extension extends to network 
(dD itself. 

Remark 7.17. Corollary 17.151 is only one consequence of Proposition 17.141 
There are other, more interesting ones (some with a graphical flavor) that 
we intend to take up in another article. For example, a network {S^,^,M} 
of rank r is weakly normal if there exist r reactions {yi -^ y'j}i=i...r with the 
following property: There is a set of distinct species y^ - {si, . . . ,5^} with 
Si € supp yi, i = 1 . . .r, such that 

det[s.-{y'^-yj)\^_^^^0. (52) 

It is not difficult to see that (15 2 p has the following interpretation: Let 
{Vi ~^ y'i}i=i...r denote the set of "reactions" obtained from {yi -^ y'i\i=i...r 
by stripping away all species not in ^*. Then fl52|) is satisfied (whereupon 
the original network is weakly normal) precisely when the resulting set of 
"reaction vectors" {y[ - yi}i=i...r is linearly independent. 

Proof of Proposition \7. 14\ Suppose that the reaction vector set {y^ - yi}i=i...r 
and the set {pi}i=i...r <= 1R_|_ satisfy the conditions of the Proposition 17.141 In 
this case, it is not difficult to see that each of these sets must be linearly 
independent. In particular, the set {y'- - yi}i=i,,,r is a basis for 5", the stoi- 
chiometric subspace for the network under consideration. 
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We begin by constructing the linear transformation Tq-. S ^ S defined by 

r 

Tocr = Y, Pi-(^ {y'i - Vi)^ Vcr € 5, (53) 

whicli can be seen to be nonsingular in the following way: Suppose on the 
contrary that there is a nonzero a* e S such that T^a* = 0. Because the set 
{yl-yi}i=i...r is linearly independent, we must have 

Pi-a* = 0, i = l...r. (54) 

Because {y'- - yj}j=i...r is a basis for S and a* is a nonzero member of S, there 
must be ^j, j = 1 . . .r, not all zero, such that 

^* = t ^jiy'j-yj)- (55) 

Insertion of this into (1541) results in the system of r homogeneous equations 

r 

Y. Pi-{yj-yj)^j=0, i = l...r (56) 

J=l 

that must be satisfied by the set {^j}j^ir. Since the determinant of the 
matrix fHOl) is nonzero, the only solution is ^j = 0, j = 1 . . . r, which amounts 
to a contradiction. Thus, Tq is nonsingular, and, as a result, detTo ^ 0. 

It remains to be shown that the requirements of weak normality are met 
by the network {^,^,^}. For this purpose, let ^q be the set of afore- 

y 

mentioned reactions {yi -^ yl}i=i...r- Moreover, let Py^y' e M_^ be chosen to 
satisfy 

Py^y::=ey, "iy^y'e^s^o, 
where e is a small positive number. Now let T^: S -^ S be defined by 
f,a := Y. Py^y'-(rW-y) 

ij^y'e.'M 



= Toa+ e 



E y-(^{y'-y) 

.y^y'^^'^.^O 



Note that T^\e=o = Tq. Because detTo i^ and because detT^ is continuous in 
e, it follows that det T^ i^O for e sufficiently small. Thus, for sufficiently small 
e, Tg is nonsingular,whereupon the network {^,^,^} is weakly normal. D 
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Remark 7.18. Corollary 17.151 derives from Proposition 17.141 by invoking the 
special choice Pi = Vi, i = 1 . . .r. When the resulting condition in Corollary 
17.151 is satisfied, the network at hand is not only weakly normal but also 
normal. To see this, it is enough to replace pi by i/i everywhere in the proof 
of Proposition 17.141 and then invoke Definition 17.11 with g^ = 1, Vs e =5^. 

8 A Concluding Remark 

When their hypotheses are satisfied, the central theorems of this article per- 
mit one to affirm concordance of a particular reaction network from inspec- 
tion of its Species Reaction Graph and, then, to invoke all of the powerful 
dynamical consequences that concordance implies [IE]. At the same time, 
these theorems provide delicately nuanced insight into network attributes 
that give rise to concordance. And, because their hypotheses are fairly easy 
to satisfy, the theorems tell us more — that concordance in realistic chem- 
ical reaction networks is likely to be common. Moreover, the theorems also 
serve to make connections between concordance results in [TS] and earlier, 
somewhat different SR-Graph-related results contained in [Hl2|[7t[9|[T7]. 

It should be remembered, however, that computational means for direct 
determination of whether a particular network — fully open or otherwise — 
is concordant or discordant, strongly concordant or not strongly concordant, 
are already available in user-friendly, freely-provided software [T3] prepared 
in connection with [18j. In most instances, this software or a variant of it, 
will be the tool of choice. 
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Appendices 



A Proof of Propositions 15.81 and 15.10 



A.l Some graph theoretical preUminaries: Ears 

In this sub-section we describe a small amount of graph-theoretical material, 
all of which is available in [1]. In a directed graph G let F be a proper 
subgraph of G (with edge-directions inherited from G). A directed ear of F 
in G is a directed path in G whose end vertices lie in F but whose other 
vertices do not. In Figure H] we show a hypothetical source block, and we 
consider the proper subgraph with vertices Si,S2,Ri,R2, R3 and the directed 
edges connecting them. Within the full graph, the path Si ^ R4 ^ S3 ^ R2 
is a directed ear of that subgraph. 




Figure 4: A directed ear 

Three propositions will be helpful later on, all applicable to nonsepara- 
ble strongly connected directed graphs, of which a source-block in the sign- 
causality graph is an example. The first two propositions culminate in the 
third. 

Proposition A.l. |4j Let F be a (nontrivial) strongly connected nonsepa- 
rable proper subgraph of a nonseparable strongly connected directed graph G. 
Then F has a directed ear in G. 
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Proposition A. 2. [3] Let F be a strongly connected subgraph of a directed 
graph G, and let P be directed ear of F in G. Then F u P is strongly connected. 

A directed ear decomposition of a nonseparable strongly connected di- 
rected graph G is a sequence of nonseparable strongly connected directed 
subgraphs of G, say Gq, Gi, . . . , G^, such that Gj c Gj+j, i = 0, .., A; - 1 and 
such that 

1. Go is a directed cycle, 

2. Gj+i = Gj u Pi, where Pi is a directed ear of Gj in G, z = 0, . . . , A; - 1, 

3. Gk = G. 

In effect, the following proposition tells us that a nonseparable strongly 
connected directed graph — in particular, a source block — can be built up 
by beginning with a directed cycle and successively adding directed ears. (It 
is left understood that we have in mind graphs with more than one vertex.) 

Proposition A. 3. |1] Every nonseparable strongly connected directed graph 
has a directed ear decomposition. 

A. 2 Two lemmas about R-to-R and S-to-S intersec- 
tions of even cycles 

Here we will argue that when, in the sign-causality graph or the SR Graph, 
two (not-necessarily directed) even cycles intersect in either an R-to-R path 
or an S-to-S path, then the third cycle so-formed is also even. We begin 
by considering the harder case, involving an R-to-R intersection. In the 
following lemma, when we indicate a path by a symbol such as R*AR**, the 
A does not denote a species. Rather R*AR** indicates a perhaps long path 
connecting reactions R* and R**. 

Lemma A. 4. Suppose that in the sign- causality graph (or in the SR Graph) 
there are distinct reactions R* and R** and three edge-disjoint (not-necessarily- 
directed) paths, denoted R*AR**, R*BR** , and R*GR** , connecting R* to 
R** . If the cycles R*BR**AR* and R*CR**BR* are both even, then so is 
the cycle R*CR**AR*. 
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Proof. We begin by noting that both R* and R** will each have adjacent 
to them three distinct edges. Because there are only two possible labels for 
each edge (the reactant and product complexes of the adjacent reaction), 
there must be at least one c-pair centered at R* and at least one centered at 
R**. (There might be more than one c-pair at R* or R** if all three adjacent 
edges carry the same complex label, but this can happen only in chemistries 
for which three or more species can appear in the same complex; we do not 
preclude this possibility.) 
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Figure 5: Illustration for the proof of Lemma [A. 41 
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We show the essentially distinct possibilities for distribution of c-pairs 
adjacent to R* and R**, relative to the paths R*AR**, R*BR**, and R*CR**, 
in Figure [51 (We have omitted from the figure cases in which 5*1 coincides 
with 5*2, 5*3 coincides with 5*4, or 5*5 coincides with 5*6, but analyses of these 
are not substantially different from those of the cases shown.) 

The proof essentially amounts to studying each case in Figured, beginning 
with an examination of the relative parity of the number of c-pairs within 
the species-to-species sub-paths connecting Si to 5*2, S3 to S'4, and 5*5 to 5*6. 
From these considerations, it can be established in each case that the cycle 
R*CR**AR* has an even number of c-pairs. 

Case (b) of the figure is typical: Because the cycle R*BR**AR* is even, 
because within the cycle there is a c-pair centered at R**, and because within 
the cycle there is no c-pair centered at R* , the parity of the number of c- 
pairs in the S'i-to-S'2 path must be opposite from the parity of the number of 
c-pairs within the 5'3-to-S'4 path. Because the cycle R*CR**BR* is even and 
because there are no c-pairs in the cycle centered at R* or i?**, the parity 
of the number of c-pairs in the S^-to-S^ path is identical to the parity of the 
number of c-pairs in the S'3-to-S'4 path. Thus the parity of the number of 
c-pairs in the S'i-to-S'2 path is opposite to the parity of the number of c-pairs 
in the S's-to-S'e, and we can conclude that the number of c-pairs in those 
two disjoint paths is odd. In the cycle R*CR**AR* there is precisely one 
additional c-pair — the one centered at R* — so the cycle is even. 

The remaining cases can be examined in the same way. D 

Lemma A. 5. Suppose that in the sign- causality graph (or in the SR Graph) 
there are distinct species S* and S** and three edge-disjoint (not-necessarily- 
directed) paths, denoted S*AS**, S*BS**, and S*CS**, connecting S* to 
S** . If the cycles S*BS**AS* and S*CS**BS* are both even, then so is the 
cycle S*CS** AS*. 

Proof. Proof of this lemma is similar to the proof of Lemma IA.4| but it is 
substantially simpler: Because c-pairs are centered at reaction vertices, all 
c-pairs will be interior to the paths 5**^45'**, S*BS**, and S*CS**. In fact, 
the hypothesis ensures that the number of c-pairs along all three paths will 
be even (respectively, odd) if that number is even (respectively, odd) along 
any one of them. From this it follows that the cycle S*CS**AS* is even. D 
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A. 3 Proof of Proposition 15.81 



In this subsection we provide proof of Proposition I5.8[ which is repeated 
below: 

Proposition 15.81 Suppose that, for the reaction network under consider- 
ation, the Species- Reaction Graph satisfies condition (ii) of Theorem \2.1[ 
Then, in any source-block of the sign- causality graph, at most one of the 
following can obtain: 

(i) There is a reaction vertex having more than two adjacent species vertices. 

(ii) There is a species vertex having more than two adjacent reaction vertices. 

Throughout this subsection we assume that the Species-Reaction Graph 
of the network under consideration satisfies condition (ii) of Theorem 12.11 
Because the putative source under consideration, viewed in the SR Graph, 
is a critical subgraph (Remark I5.5p . it follows that no two even cycles in the 
source have an S-to-R intersection. 

Consider a source-block G, which by its very nature is strongly connected 
and nonseparable, and let Go, Gi, . . . ,Gk he a. directed ear decomposition of 
the block ( §A.ip . If the block consists of just a single directed cycle (i.e., if 
G = Go); then each reaction is adjacent to precisely two species, and each 
species is adjacent to precisely two reactions. Thus the block is both an 
S-block and an R-block. 

Suppose, then, that the block consists of more than a single directed 
cycle. Then Gi is a directed cycle Go together with an adjoined directed ear 
Pq of that cycle residing within the source-block. If that ear joins a reaction 
vertex to a species vertex of the cycle, then it is easy to see that Gi will have 
two directed (and hence even) cycles with an S-to-R intersection, in which 
case condition (ii) of Theorem 12.11 would be violated. 

Thus, there are two possibilities: Either the ear Pq is a directed path 
joining two reactions of the original directed cycle Gq or else the ear Pq is 
a directed path that joins two species of that cycle. We shall suppose the 
former case and then argue in Lemma IA.6I that all subsequent ears in the 
ear-decomposition must also be reaction-to-reaction paths. Thus, at each 
stage of the decomposition, new edges are adjoined to old vertices only at 
R-vertices. In this case, the block can contain no species vertex adjacent to 
more than two edges. That is, the block is an S-block. 
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//, on the other hand, the ear Pq is a species-to-species path, then the 
block is an R-block by an almost identical argument, via Lemma [A.7I below, 
with the roles of species and reactions reversed.) 

Lemma A. 6. Suppose that, in a directed ear decomposition Go,Gi, . . . ,Gk 
of a source block G, no species vertex of Gi, < i < k-1, is adjacent to more 
than two edges. Then the directed ear P^, < i < k-1, joins a reaction vertex 
of Gi to another reaction vertex of Gi . 

Lemma A. 7. Suppose that, in a directed ear decomposition Go,Gi, . . . ,Gk 
of a source block G, no reaction vertex of Gi, < i < k - 1, is adjacent to 
more than two edges. Then the directed ear Pi, < i < k - 1, joins a species 
vertex of Gi to another species vertex of Gi. 

We provide only the proof of Lemma IA.6t for the proof of Lemma I A. 71 is 
the same, apart from a reversal of the roles of species and reactions. 

Proof of Lemma \A.6[ The ear Pi can connect (1) a reaction vertex of Gi with 
a species vertex of Gi, (2) a species vertex of Gi with a species vertex of Gi, 
or (3) a reaction vertex of Gi with a reaction vertex of Gi. We will eliminate 
the first two possibilities. 

(Case 1) Suppose that the directed ear Pi connects a reaction vertex of Gi, 
say R* , with a species vertex of Gi, say S* . We will suppose also that P, 
is directed from S* toward R*; if Pi is oppositely directed the argument in 
similar. See Figure El Because Gi is strongly connected there is a directed 
path (labeled A in the figure) entirely within Gi connecting S* to R* and 
also a directed path (labeled B) entirely within Gi connecting R* to S* . 

If, apart from S* and R* , these paths have no vertex in common, then 
they form a directed (and therefore even) cycle within Gi, and one of them 
(path A), when taken with the directed ear Pi, forms a second directed (and 
therefore even) cycle that has an S-to-R intersection with the first cycle. This 
constitutes a violation of condition (ii) of Theorem 12.11 See Figure |6t^a). 

Suppose, then, that paths A and B have internal vertices in common. 
See Figure [n](b)0 Adjacent to S* there is an outgoing edge on path A 
and an incoming edge on path B. Thus, the reaction vertices adjacent to 
5** on paths A and B are distinct. Now trace path B backward from S* 



^^To simplify the drawing in both Figures |6] and [7] we do not show paths A and B having 
common edges, but we do not preclude this possibility. The argument for Lemma fA.GI does 
not require that paths A and B be edge-disjoint. 
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Figure 6: Case 1 in the proof of Lemma [A. 61 

until its first intersection with path A at a vertex common to both paths. 
Because at least three edges within Gi meet in this common vertex, it must 
be a reaction vertex, which we call i?*. We denote by R^BS* the directed 
sub-path along B connecting i?# to S* and by S*AR* the directed sub-path 
along A connecting S* to i?*. By virtue of their construction, these two sub- 
paths have no internal vertices in common, so together they form a directed 
(and therefore even) cycle S*AW^BS*. Note that this cycle has the S-to-R 
intersection S*AW^ with the directed (and therefore even) cycle formed from 
path A and the ear Pj. This contradicts condition (ii) of Theorem 12.11 Thus, 
the directed ear Pi cannot connect a reaction vertex of Gi with a species 
vertex of Gi. 

(Case 2) Suppose that the directed ear Pi connects a species vertex of Gj, 
say 5"*, with a different species vertex of Gj, say S** , and that these vertices 
have been labeled so that the ear Pi is directed from S* to S** . Because 
Gi is strongly connected, there is a directed path (labeled A) within Gi 
that connects S** to S* and also a directed path (labeled B) within Gi that 
connects S* to S** . Should these two paths have internal vertices in common 
(Figure [Tl^a)) then, by an argument almost identical to one in Case 1, there 
would be a contradiction of condition (ii) of Theorem 12.11 

If, however, paths A and B have no internal vertices in common (Figure 
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Figure 7: Case 2 in the proof of Lemma IA.6I 

[Tl^b)), then this by itself does not constitute a contradiction of condition 
(ii) of Theorem 12.11 On the other hand, Gi is a nonseparable and strongly 
connected directed graph that is not merely a directed cycle, so the directed 
cycle S**AS*BS** must have a directed ear within Gj, which we label C 
(Proposition lA. 1 1 ) . Moreover, that ear can only connect two reaction vertices 
of the cycle because, by hypotheses, within Gi it is only reaction vertices that 
can be adjacent to more than two edges. (See Figure [Tl^c).) Thus the cycle 
S**AS*BS** has two disjoint directed ears, one (Fj) connecting a species 
vertex to a species vertex, and another (C) connecting a reaction vertex to 
a reaction vertex. 

Note that in Figure[71^c) we show only one of several arrangements whereby 
the directed R-to-R and S-to-S ears can join to the directed cycle. We show 
in Figure |8] the essentially different arrangements. (In Figure [H Rq and So 
are the cycle- vertices from which the ears are outgoing from the cycle, while 
Ri and Sj are the cycle- vertices to which the ears are incoming.) In each 
case there are two cycles — one containing the R-to-R ear and the other 
containing the S-to-S ear — that have an S-to-R intersection. (For readers 
with access to color these cycles are colored red and green.) Moreover, each 
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of these cycles is even, either because the cycle is directed or as a consequence 

of Lemma IA.4I or IA.5I Thus, condition (ii) of Theorem 12.11 is again violated. 

We have shown that Cases 1 and 2 cannot obtain, so Pi must in fact join 

two R- vertices of G,-. D 
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Figure 8: A directed cycle with disjoint directed R-to-R and S-to-S ears 



Taken together. Lemmas IA.6I and I A. 71 ensure that each source-block is 
either an S-block, an R-block, or a cycle, in which case it is both an S-block 
and an R-block. 
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A. 4 Proof of Proposition 15.101 

Here we prove Proposition I5.10[ which is repeated below: 

Proposition 15.101 Every (not necessarily directed) cycle that lies within an 
S-block or an R-block of a sign- causality graph source is even. 

We consider a sign-causahty graph source-block, G, which we will suppose 
to be an S-block. (Apart from a reversal of the roles of species and reactions, 
there is no difference in the proof if the source-block is an R-blockO) If G 
consists of a single cycle, that cycle is directed and therefore even. Suppose, 
then, that G consists of more than a single cycle. 

Let Go,Gi, . . . ,Gk be a directed ear decomposition of G, so that Go is a 
directed (and therefore even) cycle and Gk = G. Proceeding inductively, we 
suppose that for a given value oi i, < i < k, every cycle in G, is even. We 
want to show that Gj+i has this same property. Note that Gj+i results from 
the addition to Gi of a directed ear Pi, which must join a reaction vertex of 
Gi, say R*, to another reaction vertex of Gi, say R**. (Recall that G is an 
S-block, so a species vertex can be adjacent to no more than two edges.) 

Any cycle in Gj+i that is entirely within Gi is already even by the inductive 
hypothesis. A cycle in Gj+i that is not entirely within Gi must contain the ear 
Pi and is the union of Pi with a not-necessarily-directed path in Gi, denoted 
A, that connects R* to R** and that is (necessarily) edge-disjoint from P^. 
Our aim is to show that the cycle ^'1' := R*PiR**AR* is even. 

Hereafter, if G is a path in G connecting R* and R**, and if R and R' 
are vertices along that path, we denote by RGR' the sub-path of G that 
connects R with R'. No direction is implied. 

Because G, is strongly connected, there is a directed path, labeled Bq, 
that lies entirely within G, and connects R** to R*. (See Figure [9t^a).) Note 
that the cycle ^o •= R* PiR** BqR* is even because it is directed. If the path 
Bq has no edge in common with the path A, then the cycle R*AR**BqR*, 
which lies entirely within Gi is also even by the inductive hypothesis. In this 
case. Lemma [A. 41 ensures that the cycle ^'^ - R* PiR** AR* is even, which is 
what we wanted to show. 

On the other hand, suppose that the path Bq does have one or more 



""^^For an R-block, however, there is a substantiahy simpler argument. In an R-block 
every reaction vertex is adjacent to precisely two species vertices. From this and the fact 
that the block is strongly connected it follows that every cycle in the block is the union of 
causal units, so each must be even by an argument given in ^ 
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edges in common with the path A. If Bq coincides with A then the cycle 
"^t = R*P-R**AR* coincides with the even cycle ^o = R* PiR** BqR* , in which 
case we are finished. 



(a) 



(c) 



R*=R', 



R** = R, 



(b) 




(d) 



Pi : G, 




Figure 9: Schematic illustration for proof of Proposition 15.101 

Suppose, therefore, that the path A contains an edge that is not contained 
in path Bq, so the cycle "^t is different from the even cycle ^o- In this case 
we will argue that there is a sequence of cycles '^0)'^i)'^27 • • • i^ with the 
following properties: 

(i) '^0 is even, 9 = 0, . . . ,p. 

(ii) ^0 is the union of Pi and a path B0 that connects R** with R* and that 
lies entirely within Gj. 

(iii) The path B0 contains a sub-path R*BqR0 that lies entirely within path 
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A] that is, R*BgRe = R*ARe^ Moreover, for ^ > 0, R*Be-iRe-i (= 
R* ARe-i) is properly contained within R*BgRQ (= R*ARg). 

(iv) Bp = A (so that "^p = ^t^ whereupon ^t is even). 

Because of the construction of ^o it satisfies conditions (i)-(iii). Suppose, 
then, that for a particular 6 > 0, "ifg satisfies conditions (i)-(iii). We want to 
show that if path Bg is not identical to path A, then there is a cycle ^g+i 
(and, in particular, a path Bg+i) such that the requirements of (i)-(iii) are 
met. 

Suppose that path Bg is not identical to path A. Moving away from R* 
along path A, choose the first edge in A that is not in Bg, and consider 
the longest sub-path of A that contains this edge and that has no edge in 
common with Bg. It is not difficult to see that the end- vertices of this sub- 
path must be members of both A and Bg and that each such end-vertex 
will be adjacent to at least three edges of Gj+i. Because G is an S-block, 
those end- vertices must be reaction vertices, denoted Rg+i and R'g_^i, with 
Rg+i denoting the end-vertex furthest from R* along path A. (See Figure 
|9t^a)-(c).) Note that the sub-path R*BgRg (= R*ARg) must be contained 
properly within R* ARg+i. Were it not, then R'^^-^ARg^i would lie entirely in 
Bg, contradicting what has been assumed. 

Because no edge of Rg^iAR'^^^ is an edge of Bg, the sub-paths Rg+iBgR'g^^ 
and Rg+iAR'g^^ are edge-disjoint. Together they form a cycle, 

^e •= Re+iBgRg^iARg^i, 

that is contained entirely within Gj and is therefore even by the inductive 
hypothesis. 

Note that the even cycle '^g is not only the union of Pi and Bg, it is 
also he union of the path R'^^^BgRg+i and its complementary path within 
^g, the one connecting Rg+i to R'g^i via path P,. This last path we denote 
Rg+iDRg^^. Thus, we can write '^g - R'^^^BgRg^iD R'^^^. Apart from vertices 
and edges that lie in Fj, all vertices and edges in Rg+iDR'g^^ reside in Bg. 

Because the paths Rg+iAR'^^^, Rg^iBgR'g^^ and Rg+iDR'^^^ are edge-disjoint 
and because the cycles '^g and ^g are even, it follows from Lemma IA.4I that 
the cycle '^e+i := R'g^^ARg+iDR'g^^ is even. Thus, the cycle ^g+i satisfies 
condition (i) above. 



""^^For = 0, the sub-path mentioned might be the trivial one containing only the vertex 
R*, in which case Ro = R* 
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Note that cycle ^g+i is union of Pi with the path 

Be+i '■= R* BqRq^^ARq^iBqR** 
= R ARqj^^ARq^idqR 
= R ARq+iBqR 

Thus, ^e+i satisfies condition (ii). Note too that, as indicated earher, R*BgRg 
(= R*ARg) is properly contained within R*ARe+i. Thus, condition (iii) is 
also satisfied. 

So long as Eg is not identical to A these iterations can continue, with 
every new Bg^i obtained by replacing a sub-path R*BgRg+i of Bg with a 
sub-path R*ARg+i of A, with each such sub-path of A strictly longer than 
its predecessor, R*ARg. Because A has a finite number of edges, the process 
must terminate with some Bp identical to A. 

Thus, the cycle "io^ which is identical to ^, is even. 



B Proof of Proposition I5.12L 



In this appendix we prove Proposition I5.12[ which is repeated below. The 
argument here has interesting connections to the mathematics in [11], a paper 
about detailed balance in mass action systems. 

Proposition [57T21 Suppose that, in a sign- causality graph source, an R-block 
has species set y*, and suppose that no directed cycle in the block is stoi- 
chiometrically expansive relative to the -^ orientation. Then there is a set of 
positive numbers {Ms}s€y* such that, for each causal unit s -^ R -^ s' in the 
block, the following relation is satisfied: 

fR^s'M,,-e,^RM,<0. (B.l) 

Proof. Hereafter we let '^ denote the set of causal units in the R-block 
under consideration. Given the hypothesis, we want to show the existence of 
positive numbers {Msj^e^* such that 

^ < ££Z£, S/s^R^s'^'^ (B.2) 

or, equivalently, 

lnM,/-lnM, < In-^^^, y s -^ R -^ s' e '^ . (B.3) 
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Now let {ijJs}s£.y* be the standard basis for M-^*, let {€s^r^s'}s'^R'^s' ^'i/ 
be the standard basis for W^ , let T : R-^* -^ W^ be defined by 

Tq:= Y, (uJs'-uJs)-qes^R^s', (B.4) 

and let z e M"^ be defined by 

z-^ Z! hi Y^^^—es^R^s'- (B.5) 

The existence of positive numbers {Ms}s^y* satisfying f IB.SP is readily seen 
to be equivalent to the existence of g* € M'^* satisfying 

Tq* < z, (B.6) 

for then we can take Ms = exp q*, Vs € ^* to meet the requirements of (IB.Sp . 
However, by a theorem of Gale [I3l p. 46], the existence of such a q* is 

equivalent to the non-existence oi p e kerT'^n M_^ satisfying 

p-;z<0. (B.7) 

In our case, T^ : M*^*^ -^ M-^* is given by 

T'^X:= Y, Xs-~,R-~,s'i0Js' -(^s)- (B.8) 

We will say that c € M^ is a directed cycle vector if, in the R-block under 
study, there is a directed cycle such that Cs'^r^s = I ii s' ^ R ^ s is a causal 
unit in the cycle and is zero otherwise. It is easy to see that every directed 

cycle vector is a member of ker T^ n R^ . In fact, as we show in Lemma [B. II 

i^ 

below, every nonzero member of ker T^ n R_|_ is a positive linear combination 

of directed cycle vectors. 

<^ 

Suppose that p e kerT^n R_^ satisfies condition (IB. 71) . and let 

k 

p=Y,<^eCe (B.9) 

61=1 

be a representation of p as a positive linear combination of directed cycle 
vectors. Then we must have 

k 

p-z = ^a0(ce-2) <0, (B.IO) 

0=1 
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where all the ag are positive. On the other hand, (IB.lOp is contradicted 
by the fact that each cg ■ z is non-negative: To see this, suppose that, for a 
particular 9, the directed cycle corresponding to cg is 

Si -^ Rl ^ S2 -^ R2--- "^ Sn ^ Rn ^ Si. (B-H) 

Then 

cg-z = f^ln ^'''^'^' - - In n •^'^''^'''' ^ 0> (B-12) 

with the last inequality coming from the fact that no directed cycle in the 
R-block is stoichiometrically expansive. Apart from proof of Lemma IB.H 
which appears below, this concludes the proof of Proposition 15.121 D 

a// 

Lemma B.l. Every nonzero member o/kerT^ n M^ is a positive linear 
combination of directed cycle vectors. 

Proof. Here we use the terminology and results in [131 Chapter 2]. A vector 

f/ 

in ker T'^ n R^ is extreme if it cannot be written as the sum of two linearly 

f/ (jj/ 

independent vectors of ker T^ n M^ . Because ker T^ n R^ is a pointed finite 

cone (sometimes called a pointed polyhedral cone), any vector within it is 

i^ 

the sum of extreme vectors of ker T'^ n R_|_ . We will be finished if we can 

^/ 

show that every nonzero extreme vector of ker T^ n R_^ is a positive multiple 

of a directed cycle vector. 

Suppose, then, that x* is a nonzero extreme vector of kerT^n R_|_ . Thus, 

X* must be non-negative and satisfy the equation 

T^x*= Y. xUji^A^s'-uJs) = 0. (B.13) 

Let ^* denote the support of x* , and let Rf^ be the linear subspace of 
R'^'^ consisting of all vectors whose support is a subset of '^*. That is, 

Mf := {x € M'^ : Xs^R^s' = when s^R^s'i ^*}. (B.14) 

Moreover, let R^^ denote the set of vectors of Rf^ having exclusively non- 
negative components. We denote by Tj : M.f -^ R-^ the linear map defined 
by 

T'^X:^ Z Xs^R^s'{(^s'-U}s)- (B.15) 
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Because the support of x* is '^*, it is easy to see that Tjx* = T'^x* = 0. 
Thus, the dimension of ker Tj is at least one. 

In fact if, as has been supposed, x* is an extreme vector of ker T^ n ]R_|_ , 

then the dimension of kerTj must be precisely one: If not, then there is 

nonzero vector d € kerTj c M.f that is not collinear with x*. By choosing 

e > sufficiently small, we can construct another vector x^ := x* + ed that 

also resides in kerTj and that has xl^j^^^, > for each s^R^s'e'^*. 

Because d is not collinear with x*, neither is x^. By choosing p> sufficiently 

small, the vector x^ := x* - px^ e ker Tj can also be made to have positive 

<2/ 

components on '^*. Since both Xi and X2 are members of M+^(and hence 

f/ 

of R_,_ ) and since both are members of kerTj (and hence of kerT^), both 

are members of kerT^ n R_^ . Moreover, since x* = pxi + X2, x* is the sum 

of two vectors in ker T^ n R^ that are not collinear with it (and, therefore, 
not collinear with each other). This implies that x* is not extreme, which 
contradicts what we had supposed. Thus, dim ker Tj = 1. 

Now Tj can be identified with the incidence map [3] of a directed graph 
G having as its vertices the species of '^* and as edges the directed causal 
units contained in '^*. It follows from standard graph theoretical arguments 
— see Chapters 4 and 5 in [3] — that the kernel of the incidence map of a 
directed graph (when it is singular) has a basis consisting of cycle vectors. 
Each vector c of such a basis has support on the edges of a distinct (not 
necessarily directed) cycle of the graph and has components related to the 
cycle in following way: The cycle is given an orientation (e.g. "clockwise") 
and the component of c corresponding to a particular edge in the cycle is 
+1 [respectively, -1] if the edge's direction agrees [disagrees] with the chosen 
orientation. 

Because dim ker Tj = 1, every vector in kerTj is a multiple of a single 
cycle vector c*, corresponding to a cycle in G made up of edges (causal 
units) in '^*. Because x* is a member kerTj, it is a non-zero multiple of c*. 
Moreover, because the components of x* corresponding to members of '^* 
are all positive, it must be the case that all components of c* corresponding to 
members of '^^ also have the same sign, which is to say that the directions 
of members of '^^ are consistent with a fixed cycle orientation (clockwise 
or counterclockwise). That is, the cycle in G formed by members of '^* is 
directed. Thus, we can write x* = ac*, where a is positive and c* is the 
directed cycle vector corresponding to the underlying cycle in G, taken with 
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its native -^ orientation. D 



C Proof of Lemma 6.2 



Here we prove Lemma [6.2[ which is repeated below: 

Lemma 16.21 Suppose that a fully open reaction network with true chemistry 
{^,^,M} is not strongly concordant. Then there is another true chemistry 
{^,'t^,^} whose fully open extension is discordant and whose SR Graph is 
identical to a subgraph of the SR Graph for {J^,'^,^} , apart perhaps from 
changes in certain arrow directions within the reaction vertices. 

Suppose that, for a fully open network, the true chemistry {y,^,I%} is 
not strongly concordant. We let ^aug denote the set of all reactions in the 
fully open network whose true chemistry is {y,^,^}, including synthesis- 
degradation reactions of the form -^ s and s ->■ 0. To say that the fully 
open network is not strongly concordant is to say that there is a non-zero 
a e R-^ and numbers {Qiy^y'}y^y'e^^^g such that 

^ ay^y,{y'-y) = 0, (C.l) 

and, in addition, 

(i) For each y ->■ y' € S^aug such that oty^yi > there exists a species s e 
supp {y - y') for which sgn (Ts - sgn {y - y')s. 

(ii) For each y ^ y' e I%aug such that cxy^y^ < there exists a species s e 
supp (y - y') for which sgn (Tg - -sgn (y - y')^. 

(iii) For each y -^ y' ^ ^aug such that a.y^yi - 0, either (a) 0-^ = for 
all s e supp 2/, or (b) there exist species s,s' e supp (y - y') for which 
sgn Gs = sgn (y - y')s and sgn Os' = -sgn {y - y')s'. 

Our aim will be to construct a true chemistry {^,^,,^} with a dis- 
cordant fully open extension and with !^ resulting from removal of certain 
reactions in ^, retention of certain other reactions in ^, and replacement 
of still other reactions in ^ by their reverses. It is easy to see that the 
resulting true chemistry {^,'^,i^} will have an SR Graph that is identical 
to a subgraph of the SR Graph for {S^ ,^,^} , up to modification of arrow 
directions within the reaction vertices. 
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To produce the discordance we will want to show that, with ^aug denoting 
the reaction set M augmented by the degradation reactions {s ^ 0, Vs € ^}, 
there exist numbers {a?;^i/'}j/^j/'e,#„„ satisfying 

Y ay^y'{y'-y) = Q. (C.2) 

y^y'^Slaug 

and, in addition (with a as before), 

(i) For each y ^ y' e ^aug such that a^^y' i' 0, suppy contains a species s 
for which sgno"s = sgnaj,^^,'. 

(ii) For each y ^ y' e ^aug such that a^^y' = 0, cr^ = for all s e supp y or 
else supp?/ contains species sand s' for which sgn cr^ = -sgn as', both 
not zero. 

The resulting fully open reaction network, with reaction set ^aug, will then 
be discordant. (Recall Definition 14. ![ keeping in mind that the stoichiometric 
subspace of the fully open network is M-^'.) 

The construction of ^a^g, which draws on conditions (i) - (iii) just below 
(IC.ip . proceeds through sequential retention, removal, or reversal of reactions 
in J^aug in a procedure we will now describe, along with an indication of how 
the numbers {Q^y^y'}y^y'^,^^^ are to be assigned: 

1. For each y -^ y' e ^ that is irreversible do the following: 

(A) If ay^y' ^ and there exists s e supp?/ with sgnag = sgn ay^y' retain 
y -^ y' and set oiy^y' - cty^y'- 

(B) If ay^yi 4^ and there does not exist s e supp y with sgn (jg = sgn ay^y' 
(in which case there exists s' e supp?/' with sgna^' = -sgnay^y^) replace 
y -> ?/' by ?/' ^ y, and set ay'^y = -ay^yf. 

(C) If ay^yi = remove y -^ y'. 

2. For each reversible reaction pair y :^ y' e ^ do the following: 

(A) If Qy^y' 4^ OLyi^y aud thc couiplexes have been labeled so that |ttjy-*jy'| > 
I'^j/'^j/l (whereupon sgn(ay^y' - ay'^y) = sgnay^y') then 
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(i) If there exists s e suppy with sgn cTs = sgn ay^y' then retain y -^ y', 
remove y' -> y, and set ay^y' = dy^yi - dy'^y. 

(ii) If there does not exist s £ supp y with sgn ag = sgn cty^y' (in which 
case there exists s' e suppy' with sgn cr^' = -sgna^.^^') then retain 
y' -y y, remove y -^ y', and set ay'^y = dy'^y - dy^y'- 

(B) If ay^yi = cty'^y then remove both y -^ y' and y' -^ y. 

3. For each degradation reaction s ^ that is irreversible in I^aug, retain 
s -> and set as^o = cts^o- 

4. For each reversible pair s :^ e ^aug, do the following: 

(A) If ag^o 4^ (in which case Oq^s = or sgn a^^g - -sgn Os^o)) retain s ^ 0, 
remove -^ s, and set Os^o = «s^o ^ «o^s- 

(B) If Os^o = (in which case ao^s = 0), retain s ^ 0, remove ^ s, and set 
a^^o = 0. 



That f lC.2p is satisfied follows from the construction and f IC.ip . That the 
two conditions just below (IC.2p are satisfied also follows from the construc- 
tion. 
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